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I. INTRODUCTION 



The loop quantum cosmology (LQC) in the spatially flat, isotropic model was recently 
investigated in detail [l|. (For a brief review, see 0].) That investigation introduced a 
conceptual framework to construct the physical sector and gave the technical tools to sys- 
tematically explore the effects of quantum geometry both on the gravitational and matter 
sectors. With a massless scalar field as the matter source, the analysis led to the signifi- 
cant results: (i) the scalar field was shown to serve as an internal clock, thereby providing 
an explicit realization of the "emergent time" idea; (ii) the physical Hilbert space, Dirac 
observables and semi-classical states were constructed rigorously; and (iii) the Hamiltonian 
constraint was solved numerically, showing that in the backward evolution of states which 
are semi-classical at late times, the big bang is replaced by a quantum bounce. Furthermore, 
thanks to the non-perturbative, background independent methods, the quantum evolution 
is deterministic across the deep Planck regime. 

These results are attractive and affirm the assertion that the quantum geometry effects of 
loop quantum gravity (LQG) hold a key to many long standing questions [l, 0, H, H, 0]- The 
detailed investigation enriches our understanding of quantum geometry effects both in semi- 
classical and Planck regimes. However, it has a serious drawback: the critical value pcrit of 
the matter density at which the bounce occurs can be made arbitrarily small by increasing 
the momentum p^ of the scalar filed (which is a constant of motion). This is seriously 
problematic since large values of p^ are permissible and even preferred for semi-classical 
states at late times and it is physically unreasonable to expect that quantum corrections 
significantly modify the classical predictions at low matter densities. 

A careful consideration urges that this unpleasant feature is just the artifact of approxi- 
mation schemes and should be overcome by modifying the definition of the constraint oper- 
ator in a more sophisticated way. In [8|, the new Hamiltonian constraint was constructed, 
based on the same principles used in [1, [2| but improved by a more direct implementation 
of the underlying physical ideas. The new method retains all the attractive features while 
resolves the major weakness: in the improved dynamics, the big bounce occurs precisely when 
the matter density enters the Planck regime, regardless of the value of p,^ and other initial 
conditions. 

The investigations [l| and ^ altogether provide a solid foundation for analyzing physical 
issues of LQC. However, one of the limitations in this formulation is that it is obtained in a 
specific mini-superspace model and thus the domain of validity is restricted. In order to test 
the robustness of this framework, it is necessary to apply the methodology to more general 
situations. The next step is to consider Bianchi I model, the simplest case of homogeneous, 
anisotropic models. In classical general relativity, the generic approach to the initial singu- 
larity is usually understood in terms of the Belinsky-Khalatnikov-Lifshitz (BKL) scenario. 
In this scenario, along with the Bianchi IX model, the Kasner solution of vacuum Bianchi I 
model plays a pivotal role for the infinite fragmentation of the homogeneous patches close 
to the singularity j^, [13]. Therefore, a detailed construction of Bianchi I model (with or 
without matter source) in quantum theory could provide a tool to explore the quantum 
effect near the classical singularity in generic cases. Moreover, a better understanding of 
anisotropic Bianchi I models may also shed light on how the isotropy appears in anisotropy. 
For example, we can ask what kinds of semi-classical (coherent) states in Bianchi I model 
give rise to the isotropic symmetry, either exactly or approximately. Since so far there is no 
systematic procedure which leads us to the dynamics of the symmetry reduced theory from 
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that of the full theory of LQG, understanding the isotropy in the anisotropic framework 
could give a useful insight into symmetry reduction. 

In anisotropic models, the main structural difference from the isotropic case is that the 
operator 9 used in [ll, 0, 0] is no longer positive definite. However, a detailed analysis 
shows that what matters is just the operator B+, obtained by projecting the action of B to 
the positive eigenspace in its spectral decomposition. Therefore, most analytical structures 
should go through without a major modification. 

The purpose of this paper is to give a detailed analytical investigation of LQC with a 
massless scalar field in Bianchi I model. Both the precursor strategy (""/^-scheme") used 
in [H and the improved strategy ( "/i-scheme" ) used in Q are apphed and reconstructed for 
the Bianchi I model. The analytical issues are investigated in detail and the parallel results 
are reproduced. In particular, the scalar field again serves as an internal clock, the physical 
Hilbert space with the inner product is constructed such that Dirac observables are all self- 
adjoint and thereby the forms of expectation values and uncertainties of Dirac observables 
are formulated. Meanwhile, the corresponding Wheeler-DeWitt (WDW) theory is studied 
and by analogy the semi-classical states for LQC are constructed in both schemes, which 
match the semi-classical states for WDW at late times. 

However, since Bianchi I models have more variety, there are new features significantly 
different from those in isotropic models and more careful consideration is needed. In the 
classical dynamics, apart from the solutions which have the similar characteristics to those 
in isotropic models (i.e., expanding or contracting in all three directions), the Bianchi I 
model with a massless scalar also admits "Kasner-like" solutions (namely, expanding in two 
directions while contracting in the other direction or the opposite way around) . The Kasner- 
like solution with two expanding and one contracting directions eventually encounters the 
"Kasner-like singularity" (a given regular cubical cell stretches as an infinitely long line) 
in the backward evolution. When the quantum correction is taken into account, according 
to the implications from the decoupling of the problematic state in the difference evolution 
equation, we anticipate that the singularity (Kasner-like or not) is resolved and might be 
replaced by the big bounce. In the Kasner-like case, the big bounce could occur as many as 
three times, not just once, whenever the corresponding classical solution undergoes vanishing 
behavior of the area scale factors. Furthermore, in /i-scheme, the bounce is expected to 
happen at the moment when the directional factor gj := p^/lpif (which plays the same role 
as the matter density does in the isotropic case) in any of the three diagonal directions [pi 
is the area scale factor in the direction /) approaches the Planck regime of 0{hipi) {ipi is 
the Planck length), irrespective of the value of p^. In order to conclusively show that this 
anticipation is indeed the case, we need to conduct a thorough numerical investigation to 
reveal more phenomenological ramifications due to quantum geometry effects in the Bianchi 
I model. With the analytical foundation at hand, we are ready to investigate the effective 



dynamics and numerics, which we wish to study in more detail in a separate paper [11 
(Some results for the effective dynamics have been studied in (l^.) 

This paper is organized as follows. In Section [TTl we summarize the results of the classical 
dynamics of the anisotropic Bianchi I cosmology with a massless scalar source. In Section UTTl 
the Bianchi models are expressed in terms of Ashtekar variables and the phase space of the 
diagonal Bianchi models is constructed. In particular, the Hamiltonian constraint in Bianchi 
I model is given in terms of the diagonal variables and pi and serves as the staring point 
of the whole formulation. Section [IV] deals with the kinematical structure, highlighting the 
construction using holonomies to replace connection variables. The gravitational part of 
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Hamiltonian constraint is derived in Section |V] with the quantum discreteness imposed by 
^i-scheme. The WDW theory is studied in Section I VI I to show that the singularity is not 
resolved with WDW equation and to offer a systematic method to construct the physical 
sector, which then can be heuristically modified for the context of LQC In Section fVIIl the 
analytical issues in LQC are studied in detail: the total Hamiltonian constraint is derived 
and gives the evolution as a difference equation; the state in the kinematical Hilbert space 
associated with the classical singularity is completely decoupled in the difference evolution 
equation; and finally the Dirac observables are recognized and the physical Hilbert space 
is constructed with the physical inner product. Section IIVI - Section IVIII complete the 
analytical investigation for the precursor strategy (°yU-scheme). The organization of this 
part deliberately parallels that of f]\ and the notations are also used in the similar manner. 
The procedures in Section |V]- Section IVlII are then repeated all over again with appropriate 
modifications for the improved strategy (/2-scheme) in Section IVIIII Finally, the results of 
this paper and outlooks for further research are discussed in Section [IXl 

Some technical details and related topics are supplemented in the appendices. The de- 
tailed calculation for the WDW semi-classical states is shown in Appendix \M The issues 
to obtain the formulation for the isotropic model by direct inspection from the anisotropic 
results are discussed in Appendix [Bl The comparison with the results of [l|, HI, @] and their 
subtle differences are also remarked. 



II. CLASSICAL DYNAMICS OF THE BIANCHI I MODEL 



The Bianchi I model without matter source (Kasner model) has been quantized in the 
context of LQC in [5|, as has the isotropic model with massless scalar source in [H, 0, [1]. In 
order to study the LQC in Bianchi I model in the presence of massless scalar field, we first 
have to figure out its classical dynamics and solutions. This has been done in [l^, and the 
results are briefly outlined as follows. 

The Hilbert action for the cosmology with a massless scalar filed minimally coupled to 
the gravity is given by 



S — s. 



grav 

1 

8^ 



(2.1) 



In the Bianchi I model, we can take the spacetime metric to be the form 
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and the homogeneous field 0(x, t) = (pit) to be independent of the spatial coordinates. We 
then have the Hamiltonian constraint 
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where tTj are the conjugate momenta of and is the conjugate momentum of 0, which 
satisfy the canonical relations 

{ai,Tij} = 5ij and = 1. (2.4) 

The general solution is given by 

am = «... (f ) " = a.. + ''fat'^.^CO+r; | _ ^^^^ 



and 



V AtiG 



where r is the proper time defined through dr = n{t)dt, X := n = ^ is a positive but 
otherwise arbitrary function of t and the constants Hi, satisfy the constraint: 

Ki + K2 + = 1 and + ^2 + Kg + = 1. (2.7) 

Since the Hamiltonian is independent of 0, the momentum is a constant of motion, which 
is given by 

P<l> = . ;= hVSnGICK^. (2.8) 

To V oTTG 

The complete set of solutions will be over-counted unless we fix the constant ratio 
0'i,oO'2fiO'3,o/'To] thereby we introduce the dimensionless constant /C, which is to be fixed. ^ 
Note that is a monotonic function of t and thus it can be used as the internal time. 
Expressed in terms of the internal time 0, the solution reads as 

a.(0) = a.o e^^^^-^°\ (2.9) 

where the dependence on the coordinate time t is completely eliminated and only the corre- 
lation between the dynamical variables and (p is truly physical! This realizes the idea of 
emergent time [13] and shows that one can employ some subset or combination of dynamical 
fields as an internal clock [l^ . [Later in Section IVI At Section IVIIBI and Section IVIII Dt 
it will be shown that is well-suited to be emergent time also in the quantum theory (for 
both WDW theory and LQC).] 



Remark: If without any matter sources, the classical dynamics of vacuum Bianchi I model yields 
the standard Kasner solution: 

ds^ = -dt^ + t^'^^dxl + t^'^^dxl + t^^'dxl (2.10) 

where ki,K2,K3 are constants subject to the constraint: 

Ki + K2 + 1^3 = ^ and + k| + K3 = 1. (2-11) 



The choice of K, is in fact a gauge fixing. Change of JC is equivalent to the rescahng: Ki JCni and 
K-Kcf, with the new constraint: ki + K2 + = K. and = K? . The rescaled Ki and 

still yield the same evolution equation with respect to the internal time 4> as the rescaling does not 
change (1^ . 
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With the exception of the trivial Minkowskian solutions (where two of Ki vanish), all solutions 
have two of positive and the other negative, which therefore lead to the "Kasner singularity"^ 
at i = and the "planar collapse" ^at t = oo. Now with the massless scalar source introduced, the 
solution (j2.9p with the constraint (j2.7p are the modifications of (j2.10p and (j2.1ip and it also admits 
"Kasner-like" solutions (which thereby give the "Kasner-like" singularity and planar collapse) 
except that the role of t is replaced by the internal time (p. On the other hand, because of ac^, 
we also have nontrivial solutions with Kj all positive, which are not allowed for standard Kasner 
solutions. These "Kasner-unlike" solutions have only a (Kasner-unlike) singularity'^ but no planar 
collapse and behave quite the same as the isotropic solution. In particular, the constraint with 
Ki = K2 = K3 = 1/3 reproduces the isotropic solution. □ 



III. BIANCHI MODELS IN ASHTEKAR VARIABLES 

In the construction of LQC, so far, quantum gravity has mostly been analyzed in mini- 
or midi-superspace models which are obtained by ignoring an infinite number of degrees of 
freedom of the full theory. That is, the canonical dynamical variables are enormously reduced 
by the symmetry and only a finite number of relevant degrees of freedom are preserved and 
quantized. This strategy has been used in [2] and formulated rigorously in [3] for the isotropic 
models; it was also extended to diagonal Bianchi class A models (the subclass of Bianchi 
types I, II, VIII and IX, or referred as "type D") in (sj. In this section, we follow and outline 



the ideas of [16|] to describe the Bianchi models in Ashtekar variables and end up with the 



Hamiltonian constraint for the Bianchi I model in terms of the diagonal variables and pi. 



A. The Description with Ashtekar Variables 

Bianchi cosmologies are homogeneous cosmological models, in which there is a foliation 
of spacetime: M = E x i? such that S is space-like and there is a transitive isometry group 
freely acting on S. This means that we have the Killing fields on E. We can choose a 
fiducial triad of vectors and a fiducial co-triad of covectors l^;^ that are left invariant by 
the action of the Killing fields. The fiducial 3-metric is given by the co-triad "tu^: 

V = K°^^%- (3.1) 
These bases and tul satisfy the relations: 

= a/ °e^, (3.2) 

2D[M] = -Cjk'''u;rujl (3.3) 

^ At the "Kasner(-like) singularity" , a given regular cubical cell is stretched as a linear line, infinitely long 
in one direction and infinitely thin in the other two. 

^ "Planar collapse" means that a given regular cubical cell is stretched as a planar piece, infinitely large in 
two directions and infinitely thin in the other one. This is however not a spacetime singularity, since the 
spacetime curvature goes to zero (instead of infinity) when Kanser-like solutions evolve to the far future. 

^ At the "Kasner-unlike singularity", a given regular cubical cell contracts to a point. 
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where Cij^ are the structure constants of the above mentioned isometry group. The classifi- 
cation of these structure constants leads to a classification of the homogeneous spaces. The 
case that we study in this paper is Bianchi I type, in which the structure constants Cjj^ all 
vanish. 

The Ashtekar variables consist of the canonical pairs E i and Aa . The variables E j are 
densitized triads and the connections Aa are defined by 

A: = V\-^Ka\ (3.4) 

where V\ = e^j^^F^^ with F;^^ = jel{dael—dbe{+e^^ eaidbe[) the spin connection of the physical 
triad e°; Ka = KabE'^^ = q^^^'^KabE''^ with Kab the extrinsic curvature of the 3-surface and 
q the determinant of the 3-metric; and 7 is the Barbero-Immirzi parameter. In the context 
of Bianchi models, it is more convenient to consider the variables: 

E^, := Y'^'^E",, 

A,^ := AJ%1. (3.5) 

They are essentially the same as the above canonical variables, but stripped of all the 
complicated but fixed spatial dependence enforced by the Bianchi symmetry. Conversely, 
we have 

E\ = V^E^A'], 

Aa' = A/yi (3.6) 

These matrices E^j and Ai-' depend only on time (due to spatial homogeneity) and are 
our reduced canonical variables. Thus as expected, the number of degrees of freedom (in 
the configuration space) has been reduced from 9 per space point to just 9 globally. The 
Gauss, vector and Hamiltonian constraints in terms of these reduced variables are expressed 
respectively as 

g' = C.k'E^' + G e'^^A^jE^k = 0, (3.7) 
Vk = -EUAjC.k'^ + Ge'^^EMjmAkn = 0, (3.8) 
S = e'^''Cmn''E"',E^,Apk + GE^,E^,{AjAj - A^'AJ) = 0, (3.9) 

where 7 = z is used in the last line for the Hamiltonian constraint. 



B. Phase Space in the Diagonal Bianchi Models 

For the diagonal Bianchi class A models, one can fix the gauge of the 50(3) rotation so 
that Ai' and E^j are diagonalized: 

= k\i^A,^k,^'\ Pi = APe^, A\i), (3.10) 

where AV is the 5*0(3) rotation matrix and A/ is its inverse.^ In this gauge, (13.71) and 
(13. 8p are identically satisfied and the phase space is reduced to 6-dimensional and labeled by 



^ In this paper, the Einstein convention is adopted: the indices repeated on both upper and lower scripts 
are summed. When Einstein convention is not wanted, those indices which are repeated on both scripts 
but not to be summed are parenthesized. 
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(the diagonalized entries of Ai') and pi (the diagonahzed entries of E'^j) with the single 
constraint 5 = 0. In the diagonal components the symplectic structure is given by 

{c',pj} = M^GVj'j, (3.11) 

where Vq = °Li°L2°L^ is the volume of the "unit cell" V adapted to the fiducial triad and °L/ 
are the lengths of the edges of this cell with respect to °qab- (For simplicity, we assume that 
this cell is rectangular with respect to °qah-) 

In the passage from the full to the reduced theory, we introduced a fiducial metric °qah- 
There is a freedom in rescaling the diagonal entries of this metric by constants: "qa i—^ kf "qa 
for i = 1,2, 3. Under this rescaling, the variables c^, pj transform correspondingly via 

^— s> k^j^c^^^ and pi ^-^ fc^^/E^^pi and so on. Since rescaling of the fiducial metric does 

not change physics, and pi do not have direct physical meaning and it is convenient to 
introduce the new variables 

c^ = °L(,)£(^) and pi = %2"L^Pi, P2 = %°I^3P2, p^ = %i%2pz, (3.12) 

which is independent of the choice of the fiducial metric °qab- This is equivalent to the 
relation 

In terms of the new variables, the canonical relation is 

{c',pj}=^7^^G5'j. (3.14) 



^o. ~b 



The relation between the densitized triad and the 3-metric is given by qq = S^^E iE j, 
which leads to 

Pi = (02^3) sgn(aia2a3), p2 = (aiOs) sgn(aia2a3), ps = (0102) sgn(aia2a3), (3.15) 

where are the scale factors in the metric of (12. 2p . Comparing the canonical relations (12. 4p 
and ( I3.14P with the help of (13.15^ . we have 

= A'K'yG f -^TTi — —712 — -7^3 ] sgn(aia2a3), 
\a2a3 ag a2 J 

(? = in'jG ( -— TTi + -^7r2 ^vr3 J sgn(aia2a3), 

V as aifls ai J 

= Att'jG ( — -Til — -712 + — ^TTs ) sgn(aia2a3). (3.16) 

V 0-2 «-! 0102 / 

C. The Hamiltonian Constraint in Bianchi I Models 

The gravitational part of the Hamiltonian constraint is given by 

Cgrav = = ((% - KabK'^'' + (K^f) . (3.17) 

For Bianchi I models, the spatial slice is fiat and thus the intrinsic curvature ^^^R = and 
the spin connection of the triad = 0, which follows Aa = -'jKa. Consequently, we have 

Ka' = KjE\ = ^—A,'e\ = --.I^A,'E\%i%\ (3.18) 
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and 

= ^2^2 {E"'iE''jAJA^^ - E^,E^,A^'A-i) . (3.19) 
In the gauge that A/ and E^j are diagonal, the Hamiltonian reads as 

Cgrav = {c^PlC^p2 + C^PlC^Ps + C^P2C^P^) 

2 

(cVic^P2 + cVic^Pa + c^P2C%) . (3.20) 



[Note that Cgrav in fl3.20p is proportional to S in (13. 9p up to the factor q with Cij^ = 0.] 
By the identities (I3.15p . (13.16P and = |aia2a3|, (13.20p becomes 

n ifi 2^2/^27ri7r2 27ri7r3 27r27r3 OiTT? a27r| aavrfA 

C^grav = 167r G ^ \ sgn(aia2a3), (3.21) 

\ as 02 cti 020.3 O1O3 O1O2 / 

which is exactly the same as the gravitational part in (12.31) up to an overall constant factor. 
Fixing the difference of the overall factors for Cgrav used in (12.31) and in (13.211) and referring 
to (I3.20p . we end up with the Hamiltonian constraint in terms of the diagonal variables 
and Pi: 

H = Cgrav + C(j, 

= -\ {c^piC^P2 + cViC^Ps + cV2C^P3) „ ^ , + SttG , ^ vl (3.22) 

In the new variables, the classical solution (12.90 now reads as 

PiW) = Pi,o e ^ ^ =Pi,oe ^ ' 

,/^[^i±^){4>-^,) V8;;G(i^)(<^-0o) 

P3(0) = Pi,o e ^ ^ =P3,oe ^ "-f ^ . (3.23) 



Remark: The Hamiltonian constraint (j3.22p reduces to that of the isotropic model (as given in 
[sl) if we further impose the isotropy condition = = = c and pi = P2 = P3 = P with the 
canonical structure {c,p} = 8tt^G/3. □ 



IV. KINEMATICS: QUANTIZATION 

In this section, we follow the lines in pt] and generalize the procedures to seek the quantum 
kinematics for Bianchi I models. Elementary variables and the representation of their algebra 
are obtained in Section IIV Al and Section IIVBI In Section IIV C\ the triad operators are 
constructed in terms of holonomies; the fact that the eigenvalues of the triad operators are 
bounded above is discussed, suggesting that the classical singularity could be resolved by 
the quantum geometry. 
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A. Elementary Variables 



Before quantization, we should single out the "elementary functions" on the classical 
phase space which are to be lifted to their quantum counterparts. In the full theory, the 
configuration variables are constructed from holonomies he{A) associated with edges e and 
momentum variables E{S, /), which are momenta E smeared with test fields / on 2-surfaces. 
In Bianchi I models, however, because of homogeneity, we do not need all edges e and surfaces 
S. Symmetric connection A in As can be recovered by knowing holonomies hj along the 
edges e/, which lie along straight lines in the "diagonal" directions (for / = 1,2,3) in S. 
Similarly, it is sufficient to smear triads only on the rectangular surfaces of the elementary 
cell V (with respect to °qab)- 

The SU{2) holonomy along e/ in the diagonal direction / is given by: 

hi{A,ni) := Vexp / A = exp(-^/i(7)c(^V(/)) 

Jei ^ 

= cos(^^)+2s,n(^^)r,„. (4.1) 

where 2ir/ = cr/ with oi the Pauli matrices and [ii G (— cc, oo) (and [ifLi is the oriented 
length of the edge with respect to °qah)- According to Peter- Weyl theorem, a basis on the 
Hilbert space of L2 functions on SU{2) is given by the matrix elements of the irreducible 
representations of the group. Therefore, the cylindrical functions of the symmetric connec- 
tions A are obtained by the sums of the matrix elements of the irreducible representations 
of the holonomies hj, that is 

^(^) = E E#«^^'^^"/5^(^i(^'/^i))^^''^"/^^(^2(A/^2))^^^^^ (4.2) 

As mentioned in the beginning of Section UlI Bl in terms of the diagonal variables, the vector 
constraint (13. 8p is identically satisfied and thus the symmetric theory is invariant under the 
SU{2) gauge transformation. By (14. ip . on the other hand, the holonomy hj{A,fij) is given 
by exp(— |/i(/)C^^V(/)), which is not invariant if we perform a global SU{2) rotation by o"/ 1— *• 
a'j = aj. Therefore, in order to have the cylindrical functions 5't/(2)-invariant, the only 
possibility is that all the three irreducible representations in (14. 2 p must be trivial (i.e. ji = 
32 = js = 0) and in the trivial representation, R^^^^\hi{A, hj)) is simply exp(— |/i(7)c'^'^)). 

Consequently, the algebra generated by sums of products of the matrix elements of these 
holonomies is just the algebra of almost periodic functions of . A typical element of the 
almost periodic function is written as: 

gic\c',c') = J2Qe'^'"''\ (4-3) 
j 

where j runs over a finite number of integers, fiji G M and G C. As in the terminology 
used in the full theory, we can regard a finite number of edges labeled by (j, /) as a "graph" 
and the function g{ci) as a cylindrical function with respect to that graph. Analogous to 
the space Cyl of cylindrical functions on A in the full theory, the vector space of the almost 
periodic functions is called the cylindrical functions of symmetric connections and denoted 
by Cyl^. 
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Similarly, the momentum functions E{S, f) in the full theory is replaced by E{Si, fi) in 
the homogeneous case, which are obtained by smearing E j with a constant-valued function 
// on the rectangular face Sj of the elementary cell V normal to the diagonal direction / 
(with respect to °qab)- That is 

E{Sj, fj) = [ Eg^ f^j^dx^dx' = Pi"-^As„f„ (4.4) 

J Si 

where S^^ = eabcE^i^^i = ^abc\/%Pi°^'i and Asjjj is the area of Si measured by °qab times 
the orientation factor (which depends on //). In terms of classical geometry, pi is related to 
the "physical" length of the edges and the volume of the elementary cell V via 

Li = ^API\ and V = ^/imhPsl- (4-5) 

According to fl3.14p . the only non- vanishing Poisson brackets among the elementary func- 
tions are: 

{^?(c^c^c3),J9,} =47r7G^(z^,,0)e>^-'=". (4.6) 

j 

The right hand side is again in Cyl^ and thus the space of the elementary variables is closed 
under the Poisson bracket. Note that, contrary to the full theory, the smeared momenta 
E{Si, fj) now commute with one another since they are all proportional to pi due to homo- 
geneity. Therefore, the triad representation in Bianchi I models does exist and this fact will 
be exploited in Section IIV CI 



B. Representation of the Algebra of Elementary Variables 

In the previous subsection, we have seen that the algebra of the elementary variables in 
Bianchi I models is in a perfect parallel to the case of the isotropic model studied in j^. 
The elementary variables in different diagonal directions are independent of one another 
and each direction behaves as a copy of the isotropic case. Consequently, the result in [g] 
can be generalized to the Bianchi I model in the obvious way, when we construct quantum 
kinematics and seek a representation of the algebra of the elementary variables. 

Therefore, the gravitational part of the Hilbert space for Bianchi I models is T^g^av = 
L^(]R|qj^j,, d^fiBohr), the concrete construction of which is the Cauchy completion of the space 
Cyl^ of almost periodic functions of with respect to the inner product: 

(eiP^e^^n = SfirP.- (4.7) 

(Here, the right hand side is the Kronecker delta, not the Dirac distribution.) Thus, the 
almost periodic functions A/'^(c) = Af^-^^^^^^g{c^ , , c^) := e*'^'^/^ constitute an orthonormal 
basis in Tigrav- For 9i aiid g2 in Cyl^, we have 

{9i92m = 9i{^92{^ (4.8) 
and by fl3.14p the momentum operator is represented via 

d 

= -Sni'jil^— and thus (]5/5()(c) = 47r74i ^(O/^i/)-^/^, ' (4-9) 

j 
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where the Planck length square is defined as ip^ = Gh. The basis vectors N'ft are normalized 
eigenstates of pj. Writing Mpiic) = {c\fl), we have 

pilfl) = = p^^lfl). (4.10) 

Meanwhile, f l4.5l) also gives us the operators for the physical lengths of the edges and the 
volume of the cell V: 

Lj\fl) = (47r7|/i,|)^£pi|/i) = Vl/i), (4.11) 
V\fl) = (47r7)i V|/ii^2/i3| i'klfl) = V0). (4.12) 



C. Inverse Triad Operators 



In the case that the Hamiltonian has a matter sector minimally coupled to the gravity, 
when the quantum geometry is taken into account, the inverse triad operators play a key 
role 0, 17|. For Bianchi I models, the inverse triad operators can be obtained by directly 
generalizing the result of the isotropic case in [6|. The triad coefficients sgn(p/)|p/|~^''^ 
can be expressed as the Poisson brackets {c^^-*, L(/)}, which are then to be replaced by the 
commutator (divided by ih) in quantum theory. However 



as the operator corresponding 
to the connection does not exist in the full LQG [18|, [191, [20 , there is no operator on 
^grav corresponding to and we have to further re-express the Possion bracket in terms of 
holonomies a la Thiemann's trick. 



On the reduced phase space, the triad coefficients can be written as: 



sgn(pi 



1 



tr(r,/.r-){/.r-)-\L,}) 



(4.13) 



where 



V exp 







Aj Ti dx" ] = exp (°/i(j)c(^V(7)) 



cos I ) + 2 sin ( ) Ti 



(4.14) 



is the holonomy of the connection evaluated along a line of length °fii°Lj (with respect 
to "qab) parallel to the edge of the elementary cell V in the diagonal direction I. Note that 
the expression of f l4.13p is independent of the choice of °/^/. In quantum theory, the Poisson 
brackets are replaced by the commutators. This leads to the inverse triad operators: 



i sgn(pj) 



i sgn{pi 



sm 



-cos 



f ■ ( 



(4.15) 



As in the strategy used for isotropic models, to construct inverse triad operators, there are 
two ambiguities, labeled by an half integer j and a real number ^ (0 < ^ < 1) [3, |2l[. The 
general considerations in 22|, |23| urge one to set j = 1/2. For £, a general selection criterion 
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is not available and i = 1/2 and £ = 3/4 have been used in the literature. The difference on 
this choice does not affects any qualitative features and here we use £ = 1/2 since it gives a 
simpler expression.^ 

The eigenstates of the inverse triad operator are with the eigenvalues given by 



1 



1 



47r7£2jO(i^ 



L 



L 



(4.16) 



where L^^ are the eigenvalues of the length operators defined in (14. lip and the equations 



fJ'l, fJ'-s) = 1/^1 ± >l,/^2,/i3), 



COs(-°/iiC^)|/il,/i2,/i3) 



2 (1/^1 + Ai2, /is) + l/il - /i2, Ais)) 



sm 



(o°/^lC^)l/^l'/^2,/i3) = — (lyUl + °/il,Ai2,/i3) - lyUl - ''/il,/i2,Ai3)) 



2i 



(4.17) 



as well as the corresponding ones for and are used. While the choice of "fij gives no 
difference in classical level, it changes the eigenvalues of the inverse triad operators. The 
difference due to is negligible when the operators act on the states \fl) with ^ 
but becomes significant when the values of ni are comparable to In °/i-scheme, we will 
impose = k = \/3/2 as the imprint of the fundamental discreteness in the full theory 
(see Section [V Al) . while in /i-scheme, the inverse triad operators are redefined in a refined 
way such that "fij are replaced by flj, which vary adaptively on different states \fl) (see 
Section IVmB]). 

Since the inverse triad operator is diagonal with real eigenvalues in the basis of it is 
self-adjoint in Ti.^^-^. It can also be shown that the inverse triad operator commutes with 
pj despite the presence of in (14.151) . A key property of the inverse triad operators is that 
they are all bounded above. At the value /i/ = the upper bound is obtained: 



\Pi 



(27r7)-'/' £- 



PI 



(4.18) 



The fact that the eigenvalue of l/y/pj goes to zero (instead of infinity) when /i/ ^ 
commands the matter part of the Hamiltonian constraint Cmatter(Ai) to annihilate \E'(/I, 0) 
for /ii = 0, /i2 = or /is = 0. This is an important condition that accounts for the resolution 
of the singularity (see the remark in Section FVII C\i . (More physical implications from the 
finiteness of the inverse triad operator are discussed in detail in j^.) 



V. DYNAMICS: HAMILTONIAN CONSTRAINT (VSCHEME) 

To study the quantum dynamics, we need to obtain the Hamiltonian constraint in the 
quantum theory. We follow the same procedures used in [l|, ^ to construct the gravita- 
tional constraint operator in "/i-scheme. The self-adjointness of the gravitational constraint 
operator and its WDW limit are also discussed. This section focuses only on the gravita- 
tional part of the constraint operator. The total Hamiltonian including the matter sector 
and a more complete study of the analytical issues will be explored in Section IVIIi The 
improved dynamics in /i-scheme is left in Section IVIIII 



^ Note that £ = 3/4 was used in [l| while £ = 1/2 was used in 
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A. Gravitational Constraint Operator 



In the classicall theory, the Hamiltonian constraint in the Bianchi I model is given by 
fl3.22p . However, since the corresponding operator does not exist, we cannot directly use 
the Hamiltonian of this form but should recast it in terms of holonomies of the symmetric 
connections. To bring out the close similarity of the regularization procedure to that in the 
full theory, we start from the expression of the classical constraint in the full theory, which 
is 

= f d'xNe,,,F^,,e-'E^'E'\ (5.1) 



where e := y |deti?| sgn(det i?) = sgn {P1P2P3 ) ^/\piP2P3 \ and the integral is restricted to 
the elementary cell V (of volume Vo = °Li'L2°L^ with respect to °qab)- In the second line of 
(15. ip . we exploit the fact that for spatially flat models, the two terms in the first line are 
proportional to each other (because —•yKa = Aa and F^j, = e^jkAJA^^ when the spatial 
slice is fiat). Because of homogeneity, the lapse can be assumed to be constant and from 
now on we set it to one. 

To construct the Hamiltonian constraint operator, we first follow the standard strategy 
used in gauge theories to express the curvature component F^^ in terms of holonomies. Given 
a small rectangular surface ajj parallel to the face of the cell V in the I- J plane (the two 
edges of ct/j are of lengths "fifLj and °fJ.j°Lj with respect to °qab), the "afe component" of 
the curvature can be written as 

F^^r, = y "coi U r r + 7i) . (5.2) 

The holonomy around the perimeter of ajj is: 

/,(;mx,>/) = h^;^'h^;''\h^;^'^)-\h^°/'^)-\ (5.3) 

where the holonomy along the individual edges h^'^'^ is given by fl4.14p . 

Next, following the close similarity of the regularization procedure used in the full theory 

we express the triad term e^j-fce ^E E in terms of holonomies and the Poisson 



brackets of the holonomy and the volume function. In the symmetry reduced phase space, 
this is done as: 

where V = LiL2L^ = \/\piP2P3\- Note that (15.40 is exact, while (15. 2p depends on the choice 
of 7^/ ( " °/i-ambiguity" ) . 

Taking ([ED and ([5l4]) into dSH]), we have 

4 sgn(piP2P3) 



Co; 



X J2 ^"''tr V}) + 0{c' », (5.5) 



UK 
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where the term proportional to the identity inside the parenthesis in (15.20 is dropped out 
because of the trace and because the equations e"'"^ = ^J°qe^^^ and (S^x^f^q = Vo 

are used. Note that the hmit 7^/ — in (15.51) exists and it precisely equals the classical 
counterpart Cgarv in (I3.22p . However, this limit is physically fictitious and only makes 
sense mathematically. Since c is not well-defined in Tigrav? the "continuum limit" — > 
fails to exist and thus we cannot remove the regulator in the quantum geometry of the 
reduced model. (In the full theory, on the contrary, diffeomorphism invariance ensures the 
insensitivity to the regularization procedure 24, 25, 2^. The dependance on the regulator 



is traced back to the assumption of homogeneity.) 

Since there is no natural way to express the constraint exactly in terms of the elementary 
variables, the limiting procedure is essential. To faithfully convey the underlining geometry 
of the full theory, we keep "fij explicitly and treat it as the "regulator", which refiects the 
"fundamental discreteness" of the full theory as an imprint on the reduced theory. (The 
appropriate value of is to be determined soon.) The resulting regulated constraint is: 

87r734i°/ii°/i2°/i3 y i J i J k i k . j 

= 3.2 o o o 1 /^l^ ^^^y /^2C Sm — — \/cOS — - - COS — - ^ sm — - 

7r7'^£pi°/ii°yU2 yU3 \ 2 2 2 2 

+ cyclic termsj, (5.6) 

where "cyclic terms" denote the same expression as the term(s) written explicitly inside the 
curly bracket but with the cyclic index replacements (1 — 3, 2 — 1, 3 — 2) and (1 — 2, 2 — >^ 

3,3 —> 1). (The subscript "o" used for Cgrav,o is to indicate that the ordering ambiguity of 

quantization has not been taken care of and thus Cgrav,o is not self-adjoint.) The operator 

Cgriir^o acting on the eigenstates of pi gives 

C'sral,ol/^l5/^2,/i3) 



87r734j°/ii°/i2°/i3 

X ||V'/^i,/i2,M3+°M3 ~ ^Mi./^2,M3-°M3 I (^1/^1 + 2 °/^l)/^2 + 2')i2,/^3) — 1^^1—2 /i2 + 2 ')i2, /is) 
-|yUi + 2 7ii,/i2-2°/i2, Ai3) + 1/^1-2 °/il,/i2-2>2, /is)) 

+ cychc terms I , (5.7) 
where (I4.17P is used. 

Regarded as a state in the connection representation, (every element of) the holonomy 

fo. \ — 

h\ is an eigenstate of the area operator d/ = \pi\: 

a h'^°^'\c) = (47r7 h^j^'\c). (5.8) 

Following the same strategy (7i-scheme) in Jl|, B, [pi, we demand this eigenvalue to be be A = 



2v/37r74i, the area gap of the full theory |27|, l28|,l29|. This yields 7i/ = k := A/(47r74i^ 
V^/2 and fixes the " °/i-anibiguity" . 
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B. Self-Adjointness and the WDW Limit 



In order to construct the physical sector with the physical inner product (which will be 
done in Section IVIip . it is advantageous to have the Hamiltonian constraint self-adjoint. 
Note that fl5.6p is not self-adjoint. This happens simply because of the factor ordering 
ambiguities common in quantum mechanics. There are two natural ways to resolve the 
ordering problems. The first is to take the self-adjoint part of (15.61) : 



roy = - 

grav 2 



/7(°m) I /^(°M)t 

grav,o ' grav,o 



(5.9) 



The second is the nature re-ordering adopted in [ij, in which we "re-distribute" the sin(7ic) 
terms to get the following self-adjoint constraint: 

4z sgn(piP2P3) 



x<! sm{°fiic^] 



sm( )V cos( , 



cos( )V sm{ 



+ sm(7i2C ) sm( cos( ) — cos( )\/ sm( 



sin(°/i2C^) 



sin(°yUiC^) 



-|- cyclic terms I . 



(5.10) 



■ (°uY 

When acting on the eigenstates of p/, Cgrav gives 

1/^1, At2, Ai3) 



+ 



Ml 


+ 2°M1 


M2+2°M2, 


— V^M1+2°M1 


M2 +2°M2 ! 


+ 


^Ml, 


M2. 


^1 


M2. 




+ °M3 




M3 - °M3 






M3' 


f°M3 


M3 


-°M3 


Ml 


~2°M1 


M2 +2°M2. 


V^Ml -2°M1 


M2 +2°M2, 


+ 


V^Ml 


M2- 


V^Ml 


M2. 


M3 


+ °M3 




M3 - °M3 






M3 


f °M3 


M3 


-°M3 


Ml 


+ 2°Mi 


M2 -2°H2, 


— V^Ml+2°Mi 


M2 ~2°M2, 


+ 




M2i 


V^Ml 


M2. 


M3 


+ °M3 




M3 - °M3 






M3 


f °M3 


M3 


-°M3 


Ml 


-2°M1 


M2-2°M2. 


^Ml -2°M1 


M2-2°M2. 


+ 


^Ml 


M2. 


^Ml 


M2. 


M3 


+ °M3 




M3 - °M3 






M3+ M3 


M3 


-°M3 



+ cyclic terms I 
and Cgral gives 

^grav l/LH,/^2,/^3) 



|/il + 2>i,/i2 + 2 7^2, Ai3) 
|/il-2°/ii,/i2 + 2°/i2,/i3) 
|/il + 2°/ii,/i2-2°/i2,/U3) 

|/il-2>i,/i2-2>2,/i3) 
(5.11) 



1 



167r734i°/ii°/i2°/i3 



X 



+ 



Ml + 2°Mi 
M3 + °M3 


M2 


V^Ml 

M3 


+ 2°M1 
-°M3 


M2, 


+ 


^Ml, M2 + 2°M2, 
M3 + °M3 


-^1 

M3 


M2 + 
-°M3 


2°M2 


Ml -2°M1 
M3 + °M3 


M2 


V^Ml 

M3 


-2°Mi 
-°M3 


M2, 


+ 


^Ml, M2+2°M2, 
M3 + °M3 


-V^Ml 

M3 


M2 + 
-°M3 


2°M2 


Ml +2°M1 
M3 + °M3 


M2 


^Ml 

M3 


+ 2°Mi 
-°M3 


M2. 


+ 


^Ml , M2 -2°M2. 
M3 + °M3 


-^Ml 

M3 


M2 - 
-°M3 


2°M2 


Ml -2°ni 
M3 + °M3 


M2 


V^Ml 

M3 


-2°M1 
-°M3 


M2, 


+ 


^Ml. M2 -2°M2. 
M3 + °M3 


-V^Ml 

M3 


M2 - 
-°M3 


2°M2 



+ cyclic terms I . 



|/il-2 7ii,/X2 + 2>2,/i3) 
|/il+2°/ii,/i2-2°/i2,/i3) 
|;Ul-2°/ii,/i2-2°/i2,/U3) 

(5.12) 
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In this paper, we choose the second method of symmetrization, since Cgrav has certain 
qualities better than Cgrav • In particular, in the isotropic cases, Cgrav gives positive definite 
G (to be defined later), which is a pleasing feature when the physical sector is constructed. In 
the anisotropic Bianchi I model, the positive definiteness of G no longer holds; nevertheless, 
choosing Cgrav enables us to generalize the isotropic results in an obvious way. 

The Hamiltonian constraint equation differs markedly from the WDW equation of ge- 
ometrodynamics in the Planck regime because it crucially exploits the discreteness under- 
lying quantum geometry. In the continuum limit °fii —>■ (although physically fictitious), 
one expect that the LQC constraint equation will reduce to the WDW equation. According 
to the analysis in [l| (and more details in [6|), we can follow the arguments therein and 
conclude that the gravitational part of the Hamiltonian operators reduces to: 



C'Si,o^(p) ^ 1287r24|sgn(pip2)y^^^ + cyclic terms}vl/(p) 

=: C^Z';'om, (5.13) 

C[£^ip) ^ 64vr24i{^^^)^ (^^9^ + 9^^^) + terms}vl/(p-) 



C^r^Mp) - 64vr^4i|^^^^ f^V^T^ + ^VM^) + cyclic termsU(p-) 

vIpsI \9P2 opi dpi dp2j J 

=: C^T^iP)- (5.15) 

In the above approximations, ~ stands for equality modulo terms of the order 0{')3j) (where 
'^j := 4:7r'jip°fij). That is, in the limit when the characteristic scale is much greater than that 
associated with the area gap (i.e. \pi\ ^ °pi or ^ "/ij), quantum geometry effects can 
be neglected and the LQC difference operators reduce to the WDW differential operators. 
The different ordering factors give rise to different WDW limits. In particular, Cgrav and 

" (°u) 

Cgrav differ only by a factor ordering: 

(C^;^)' - Cgl) ^(p) = ^^"^ , (sgn(pip2) + sgn(pip3) + sgn(p2P3)) ^(p). (5.16) 

The difference is negligible for large (volume) scales. The operator will be used 

extensively in Section IVll 

Note : The results obtained in this section can be formally "translated" for the isotropic case by 
inspection. The translation and the related issues of isotropy reduction from anisotropic models in 
"//-scheme are remarked in Appendix IB 1[ □ 



VI. WHEELER-DEWITT (WDW) THEORY 

In this section, we make a digression to study and construct the WDW quantum theory in 
the connection dynamics. This construction will serve two purposes: first, it will introduce 
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the key notions for constructing the physical space and enable us to generalized these notions 
heuristically for the completion of the program in LQC; second, we will be able to compare 
and contrast the results of the WDW and LQC theories, thereby extracting out the essential 
ingredient given by the quantum geometry. In particular, the WDW theory does not resolve 
the classical singularity but the semi-classical (coherent) states in the WDW theory can be 
used to build up the initial states at late times in LQC theory. 



A. Emergent Time and the General Solutions 

Recall that the classical Hamiltonian constraint in terms of the variables and pi is 
given by (13.221) : that is 

Cff) Cgrav 

SttG \piP2P3\~'^^^pI = ^ i ^^^^^^^ c^c^\/\piP2 \ + cyclic terms I . (6.1) 

In WDW theory, we have the kinematic Hilbert space TiSi?^ ~ L^(R^, dp dcf)). The operators 
Pi and (p act as multiplications while and act as differential operators, represented as: 

c^-qj = iflSn^G^— and ps'^ = -ih—. (6.2) 
opi 0(p 

To write down the quantum constraint operator, we have to make a choice of the factor 
ordering. Since our purpose is to compare the WDW theory with LQC, it is most con- 
venient to use the same factor ordering that comes from the continuum limit of the LQC 
constraint. Therefore, we chose defined in (15.151) as the gravitational part of the 

WDW Hamiltonian constraint. The total WDW Hamiltonian constraint then reads as 

S = SnG imi-'l'-^^ ( ^VI^It^ + j^VMt^) + cychc terms}* 

=: - (03(/^i> /^2) + e2(/ii, /is) + ei(/X2, /is)) ^ =: -0(/2)^, (6.3) 

where we define pi =: 47r7£pj/i/ and 5(/i) := |/ii/i2/i3p^^^. The WDW equation now has the 
same form as the Klein-Gordon equation, with playing the role of time and G of the spatial 
Laplacian. Therefore, the form of the Hamiltonian constraint in the WDW theory suggests 
that can be interpreted as emergent time. Unlike the isotropic case, however, B is not 
elliptic and hence not positive definite. Fortunately, this problem can be easily overcome. 
We will see soon that the subspace of negative eigenvalues of 9 is actually non-physical and 
it turns out only 6+, the projection on the positive spectrum of 6, does matter. 

Note that if \l/(/i, 0) is a solution of (16.31) . so is n/^(/I, 0) (where Hj are the orientation 
reversal operators, which flip ni to — /i/ for 1 = 1,2, 3) since 9 commutes with H/. Because 
the standard WDW theory deals with geometry, it is completely insensitive to the orientation 
of the triad and H/ should be considered as large gauge transformations. Furthermore, 
Hj = 1 and there are only two eigenspaces of H/, one of which is the symmetric sector and 
the other anti-symmetric. Therefore, it is natural to work with the symmetric sector. Thus, 
the physical Hilbert space will consist of the suitably regular solutions ^^(/i, 0) to (16. 3p which 
are symmetric under /x/ — + — /x/ for / = 1, 2, 3. 
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To obtain the general solution of (16.30 . we first note that the operator enclosed by the 
curly bracket in (16. 3p is self-adjoint on L|(M^,(i^/i), the invariant subspace of L^(M^,(i^/i) 
under n/=i,2,3- Therefore, 6 (and Qj as well) is self-adjoint on L|(]R^, 5(/i) cPfi). It is easy 
to verify that the eigenfunctions of can be labeled by A; G M^: 

e^{fl) = ^-1^ e^^^'-^l'^ile^'^^^l^^le^'^^^^l^^l. (6.4) 

The corresponding eigenvalues are given by 

e eg(/i) = uj^{k) ej:{fl) with uj^{k) = IGnG (^kik2 + kik^ + k2k3 + , (6.5) 

(where the factor 3/16 is an artifact of the factor ordering choice). The eigenfunctions satisfy 
the orthonormality relation: 

J d% B{fl) e^ipi) e^, {pi) = 5\k, k') , (6.6) 
(where the right hand side is the standard Dirac distribution) and the completeness relation: 

d%B{fI) ej:{fl) ^(/I) = for VA?, iff ^(/2) = (6.7) 



for any ^(/2) G Ll{R^,B{fl) d^fi). 

Unlike the case of isotropic model, however, 6 is not positive definite on \l/(/i) G 
L|(M^, 5(/2) d'^yu) and could be negative. If cj^ < 0, the frequency has imaginary part 
and the corresponding plane wave solution to (16.31) blows up in the far future (0 oo) or 
far past {(j) — *• —oo); therefore, the modes associated with negative lu"^ should be discarded 
and only those with a;^ > are allowed. The relevant space is the "positive" subspace 
L|(M^, 5(/i) (Pjj; iJ^ > 0) and what matters is just the operator 6_,_, obtained by projecting 
the action of to the positive eigenspace in its spectral decomposition.^ 

With the eigenfunctions in hand, we can express the "general" solution to (16. 3p as 

^(/i,0)= / rf='A;('l+(fc)e,^(/I)e^-<^ + ^_(^)eg(/i)e-*-^) (6.9) 



^ The physical condition > is reminiscent of the constraint (|2.7p on the classical solution parameters 
Ki, K2, ^3 and which leads to = 2(kiK2+kiK3+K2K3). If we have the identification: fci_2.3 = /Cki^2.3, 
then the frequency gives 



u;2(fc) = 8^G \K?k\ + « S^G/C^^ = ^ > 0, (6.8) 

(where the approximation is valid when ^ h^jG or equivalently /Ck^ 1.). This shows that the 
physical condition > simply means is real. For a semi-classical state, the identification A:i_2,3 = 
^^^1,2,3 is sensible, since its Fourier transformed function ^'±(fc) (to be defined below) is sharply peaked 
at ^1,2,3 = A^'«i,2,3- (See Section IVI CI for more details.) 
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for some '^±{k) in L^(R^, d^k). Following the terminology used in the Klein-Gordon theory, 
if \E'±(/i;) has support in the region of positive kj, we call the solution "incoming" (or "ex- 
panding") in the direction /, while if it has support of negative kj, it is said to be "outgoing" 
(or "contracting").^ 

If \^+/_(/c) vanishes, the solution is said to be negative/positive frequency. Thus, every 
solution (16.91) admits a natural decomposition into the positive and negative frequency parts. 
The positive (respectively negative) frequency solutions satisfy the first order (in 0) equation 
which can be regarded as the square- root of (16. 3p : 

^ = ±zv^vp±, (6.10) 

where is the square-root of the positive self-adjoint operator 9_,_ defined via spec- 

tral decomposition on L'g{M.^ , B_{fl) d^fi; u"^ > 0). Therefore, a general "initial" function 
\E'±(/i, 0o) = /±(/i) G L|(]R^,_B(/i) d^fi;uj'^ > 0) at = 0o can be "evolved" to a solution of 
(EjDby: 

vl>±(/i,0) = e±^V^(^-^°)vl>^(/2,0.). (6.11) 



B. The Physical Sector 



To endow the space of the WDW solutions with a Hilbert space structure, we follow the 
same strategy of [l[ . The idea is to introduce the operators corresponding to a complete set 
of Dirac observables and select the required inner product by demanding that those Dirac 
operators all be self-adjoint. In the classical theory, a complete set is given by and 
(/2 at some particular instant (j) = (po), since the set of these uniquely specifies a classical 
solution and thus represents a point in the classical phase space. In WDW theory, given a 
(symmetric) solution \E'(/i, 0), since p^j, commutes with 9, 



P0^(/i,0) := - 




(6.12) 



is again a (symmetric) solution to the WDW equation. So, we can just retain this definition 

of from H^^^ . On the other hand, (16.111) enables us to define the Dirac observables 
For a given (symmetric) solution \E'(/i, 0), we decompose it into the positive and negative 
parts \E'-|-(/i, 0), freeze them at = (p^, multiply them by and evolve each of them by 
fl6ll|) . That is 

/;7U^(/^,0) := e^v^('^-'^°V/l^+(/^,0o) + e-^v^(*-'^°V/l^-(/^,0o). (6.13) 
The result is again a (symmetric) solution to (16. 3p . 



The name of "expanding/contracting" makes sense for semi-classical solutions, for which fc/ is supported 
mainly in the vicinity of ICkj as mentioned in the previous footnote. According to (|2.9p . if k/ is posi- 
tive/negtive (and > 0), the scale factor a/ in the direction / is expanding/contracting with respect to 
0. [Note: The nomenclature adopted here is opposite to that in due to the minus sign used in the 
equation there: = —hy^S/lGirG k* , which, in the anisotropic notations, gives rise to an extra minus 
sign for the relation between and in (|6.20p .] 
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Note that both and /i/|</,„ preserve the positive and negative frequency sectors and 
hence we have super- selection between these two sectors. In quantum theory, we can restrict 
ourselves on one of them. From now on, we will focus on the positive frequency sector and 
thus drop the suffix + for ^! . 

Every WDW solution is uniquely determined by its initial condition 'if{fI,(j)o) and the 
Dirac observables have the actions on the initial state: 



|/i/|^(/2, 0o) and ]30^(/2, 



(6.14) 



The unique (up to an overall scaling) inner product which makes these Dirac observalbes 
self- adjoint is therefore given by 



(^l|^2) 



phy 



rf^/i5(/i)^l(/i, 0o)^2(/i, 



(6.15) 



(Note that the inner product is conserved, i.e., independent of the choice of the initial instant 
0O-) As a result, the physical Hilbert space Ti.^^ is the space of wave functions '^{Jl^cj)) 
which are symmetric under — — and have finite norm by (16.151) . The representation 
of the complete Dirac observables on ?^^y^ is given by 



and p^'${il,(j)) 



(6.16) 



With the Dirac operators as well as the physical inner product in hand, we are able to 
calculate the expectation and uncertainty of the physical observables for a given physical 
state and thereby to extract the physics in WDW theory. 

The same representation of the algebra_of Dirac observables can be obtained by the more 

' 1, m, 13] , as has been done for the isotropic 



systematic group averaging method 



30, 31 



model in pLI]. There may be some complications in anisotropic cases and we defer this issue 
to future study. 



C. The Semi-Classical (Coherent) States 

With the physical Hilbert space and a complete set of Dirac observables available, we can 
introduce semi-classical (coherent) states and study their evolution. We first fix the time 
at = 00 to construct the semi-classical state which is peaked at = and = /i^^. 
To have a semi-classical state corresponding to a universe with the characteristic quantities 
much larger than Planck scales, we demand that ^}^^ ^ and p^ ^ h\/G (or equivalently 
K,K^ ^1). At the initial time = 0o, consider the coherent state: 

.{fl,<Po)= [ ^ d^A;^,(fc)e,^(/7)e^-(^^)(<^°-<^^), (6.17) 
where the Fourier amplitude '^aik) is given by 

_ \ \ (fc]^-KKi)2 (fcg-KKj)^ (fca-KKg)^ 

^^k) = ,, ^ e (^-1)^ e (^-2)^ e (^-3)^ , (6.18) 
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which is sharply (if (7^ ^ /C/tj) peaked at = /C/?. The constants /TJ^ and are related by 

^ '^(6 / <o (2cri )^ /o \ '^d) ^ 



Ai30„ = e'^-'s)^ e ^ ^ = e(2<^i' e v '-^ ^ ^ (6.19) 

which are suggestive of the classical solution (13.231) . Similarly, referring to (12. Sp . we set 

pl = hVs^K,K^. (6.20) 

The parameters ^1,2,3 and satisfy the constraint (12. 7p . 

As the action of the Dirac observables on the initial states and the physical inner product 
are well defined in (I6.14p and (I6.15p . in Appendix|Al we are able to calculate the expectations 
and dispersions of the Dirac observables. The results are: 

• The expectation value of 

• The dispersion square of /i/|</,„: 

A^(^/Uo) = ((m/Uo)^) - {i^iUo? ■= (/^/Uo) l%^,p;)phy - (^/UJ^ 

^ {e^ -l){fx%f. (6.22) 

• The expectation value of pip: 

• The dispersion square of p^: 

2p* 



,★2 



{al{K2 + K^f + al{Ki + K^f + al{Ki + Ka)^). (6.24) 



As the evolution of the coherent state is concerned, because of (16.111) . the solution 
^/i*,p* (yU, 0) defined by the initial state p,*^ ^p*^{fI,(l)o) can be obtained simply by replac- 
ing 0o with in (I6.17p ! This implies that when the state is evolved to the instant 0, it is 
again a semi-classical state peaked at 

(/..I,) -e(^e^('^)(^°-^^ ^/.4 (6.25) 
and with the uncertainty spread: 

A'iM^ie^ (6.26) 
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Comparing fl6.19p and fl6.25p . we have 



^^i<p = ^^I<t>o e ^ ) (6-27) 

which shows that the semi-classical state continues to be peaked at the trajectory exactly 
described by a classical solution as in (13.231) . On the other hand, (p^) and A^(p^) are 
independent of time. 

The fact that the semi-classical states follow the classical trajectories implies that they 
are destined for the singularity (Kasner-like or not) in the backward evolution. In this sense, 
WDW evolution does not resolve the classical singularity. 

Remark: To get the semi-classical (coherent) state in the context of LQC, we simply replace e^{P) 
in (j6.17p with e^^\j2) (the symmetrized eigenfunctions of 0, to be defined in Section [VII Bp . Since 
e^^\fl) erlfl) for 1^/1 %i, the wavefunction given by (|6.17p is used as the initial state at late 



times to conduct the numerical analysis for LQC See Section IVlI Dl for more details of LQC 



semi-classical states. □ 



VII. ANALYTICAL ISSUES IN LQC ("/i-SCHEME) 

In this section, we come back to the model in LQC and more detailed analytical issues 
are explored. We include the matter sector to the Hamiltonian and show that the scalar field 
can again be used as emergent time. The total Hamiltonian constraint gives a difference 
revolution equation and the general solutions to it are closely examined. The results enable 
us to construct the physical Hilbert and the Dirac observables. In Section IVII C| we also 
prove that the state in the kinematical Hilbert space associated with the classical singularity 
is completely decoupled from the evolution, which adds evidence to the occurrence of the 
big bounce. 



A. The Total Hamiltonian Constraint (Matter Sector Included) 



Since = SvrGpV \/\piP2P3\ in (13.220 . the corresponding operator C*^ in LQC acting on 



the states of Ti. 



total 
kin 



n 



grav 



is given by 



According to (14.160 . this gives 
where 

Bifl) 



(47r7)2£| 



-BiiI)\lI)0 87iGpj 



PI 



X 



(7.1) 



(7.2) 



(7.3) 
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The total Hamiltonian constraint Cgral + C{ 



then yields 

S7rGB{iI)pl \fl, 0) + (47r7)i4iC'Sl 1/^, <l>) = 0- 



(7.4) 



Both the classical and the WDW theories suggest that (j) can be used as emergent time; 
i.e., is to be thought of as "time" while fl as the genuine, physical degrees of freedom which 
evolve with respect to 0. To implement the total kinematical Hilbcrt space with this idea, we 
choose the standard Schrodinger representation for but the "polymer representation" for jl, 
which captures the quantum geometry effects; that is H^^^ — L'^(M.^^^j., d^/XBohr)<H)-^^(lR, d4>)- 
In this representation, P</) = f ^ and the total Hamiltonian constraint reads as 



d 



^^-*(/I,0) = 2 oj,o^,o^, t-^^^^]"' {^3^^(^^ ^(^^ + ^2 + 2>2, Pis: <t>) 

-C+-(ll) + 2°//i, //2 - 27^2, 1^3, (t>) 

+Cr{pi) - 27^1, /i2 - 27/2, /i3, 0) 

+ cyclic terms I 

= : - (03(/il,/X2) + 62(/il,/^3) + 0l(/i2,Ai3)) 

=: -G(/I)vl/(/I,0), 



(7.5) 



where 



+ 



+ 



+ 



^/\{^JLl + 2°ni)n2{pi3 + "A^s)! - VI(a*i + 2''A*i)a*2(a*3 - °A*3; 

y|^l(^2 + 2°^2)(/U3 + ')i3)l - \/|/Wl(/W2 + 2"/X2)(/i3 " 7^3)1 
VK/^l - 27fi)At2(At3 + 7*3)1 - VK/^l - 27il)/X2(A«3 - 7*3)1 

y|/Xi(/X2 - 27/2) (/^3 + 7^3)1 - \/|/^l(/^2 - 2°//2)(/^3 - 7^3)1 
VK/^l + 27Xi)Ai2(/U3 + 7^3)1 - + 27^0/^2 (/^3 - 7^3)1 

v/|/ii(/i2 - 2<'//2)(/i3 + 7*3)1 - A/|Atl(A^2 - 27/2)(/i3 " 7*3) | 



+ 



V|(/^l - 2°H^)H2{I^3 + 7^3)1 - \/|(/^l - 27^i)/X2(/^3 - 7^3; 

V|Ati(At2 + 27(2)(At3 + 7*3)1 - a/|a*i(a*2 + 27t2)(At3 - 7^3)1 



(7.6) 



and similar for C2 and in the cyclic manner. The coefficients Cj have the following 
properties: 



C3++(0, /i2, Us) = C3+ (0, H2, /is) = C'g +(0, yU2, yUa) = C3 (0, /X2, /^s), 
C3++(/.l,0,/.3) = C3+-(/Xi,0,/.3) = C3-+(/.l,0,/.3) = C3--(/.l,0,/X3), 



(7.7) 



and 



C3++(0, 0, /is) = C3+- (0, 0, /xs) = C3-+(0, 0, /xa) = (0, 0, /X3) = 0, 

C++ ill,, 112,0) = C3+-(/Xi,/X2,0) = C3-+(/Xi,/X2,0) = C3--(/Xi,/X2,0) = 0, (7.8) 
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for V /ii,/i2,At3- (C*! 2 have the similar properties in the cyclic manner.) For later conve- 
nience, we also define 

Cf^{f,j,f,K) := —^^mr'CfHfl), (7.9) 

Z fli fl2 fl3 

which is independent of /i/ and has the properties: 

= = C^+i2"fx,,-2%) = C^-i2"fx,,2%) = 0, (7.10) 

C'3±±(0,/i2), C'3±^(/ii,0)^OO, (7.11) 

and 

Cf^ifJ'i, 7^ or CX3, elsewhere, (7-12) 

for V /ii,/i2- (Ctt have the similar properties in the cyclic manner.) Also note that Qj is 
independent of /i/. It is easy to see that B and each of individually are self-adjoint on 

L^{Rl^^^,B{fI)d%Botr). 



B. Emergent Time and the General Solutions 

The resulting Hamiltonian euqation (17. 5p is cast in the way very similar to that of the 
WDW constraint ( 16. 3p . The key difference is that the ^-independent operator O is now 
a difference operator rather than a differential operator. In the same spirit as for WDW 
constraint, the LQC quantum Hamiltonian constraint can be regarded as an "evolution 
equation" which evolves the quantum state in the emergent time 0. 

As in the isotropic model, since 9 is a difference operator, the space of physical states 
(i.e. appropriate solutions to the constraint equation) is naturally divided into sectors, each 
of which is preserved by the "evolution" and by the action of the Dirac observables. Thus, 
there is super- selection among these sectors. Let £|e^| (|e/| < °/i/) be the "lattice" of points 
{\ei\ + 2n°fij] n E Z} on the /i/-axis, be the "lattice" of points { — le/l -|- 2n°iii]n e Z} 

and C^j = C\^j\ [JC^\^j\. Also let 'H^^^_^^ _|_|^^| ^^^^^ (8 of them) and Hf^^ denote the physical 

subspaces of L^(]R|qj^j,, _B(/i) ci^/XBohr) with states whose support is restricted to the lattices 
>C±|ei| ®'^^±|e2| ®'^±|e3| a^d ® ® A3 respectively. Each of these subspaces is mapped to 
itself by 9. Furthermore, since Cgral is self-adjoint on H^^^ = L^(]R|q[^j,, d^/XBohr)) it follows 
that 9 is self-adjoint on all these physical Hilbert subspaces. 

Note however that since C\^j\ and are mapped to each other by the operators 11/, 

none of 'H^^^^^ ^^^^^ ^^^^^ but only Tif^^ is left invariant by ni^2,3- The operators 11/ reverse 
the triad orientation and thus represent a large gauge transformation. In gauge theories, 
we have to restrict ourselves to the sector corresponding to an eigenspace of the group of 
large gauge transformations. Since there are no fermions and no parity violating processes 
in our theory, we are led to choose the symmetric sector (i.e. the subspace of "Hf^'^ in which 
^'(/I) = ^(-/ii,/i2,/i3) = ^(/^i, -yU2,/i3) = ^(/ii, At2, -Ai3)), which will be of our primary 
interest . 

To explore the properties of the operator 9, Let us first consider a generic e (i.e. none 
of e/ equals or On each of the 8 Hilbert space Ti.^^^^^ _|_|^^| _|_|^^|, we can solve for the 
eigenvalue equation Qe\{fl) = (9i + 92 + 93) ex{fl) = Xe\{fl). Since 9/ is independent of 
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FIG. 1: (a): Given the values of e\-^^ at the lattice points marked •, the value at the point marked x 
can be uniquely determined via the recursion relation (|7.14|) . (b) : If the values of 6^12 are specified 
on the boundary points (marked •), all the values at the points marked x can be obtained by 
iterating the scheme depicted in (a). 

/i/, we solve the eigenvalue problem by separating the variables; that is, let eA(/xi, /i2, /^s) = 
eAi(/Ui) ex.,{n3) and solve 

03eAi2(/^i,/^2) = Ai2eAi2(Ati,yU2), for eAi2(yUi,/i2) = eAi(Aii) eA2(Ai2), 
02eAi3(/ii,/i3) = Ai3eAi3(Aii,/i3), for eAi3(/ii,/i3) = eAi(/ii) eA3(/i3), 
0ieA23(Ai2,/i3) = A23eA23(Ai2,/i3), for eA23(/i2,/^3) = eA2(/i2)eA3(/X3), (7.13) 

with A = A12 + Ai3 + A23. Each of (17.131) can be solved recursively for any arbitrary Ajj 
provided that the appropriate boundary condition is given. For example, the detail of 
03eAi2(/^i,/^2) = Ai2eAi2(AH,/^2) reads as 

C'3^"^(Aii,/i2)eAi2(/^i + 2>i,/i2 + 2°/i2) - (53"+(/ii,/i2)eAi2(/^i - 2>i,/i2 + 2°/i2) 

-C^~{fll,fl2) eAi2(AH + 2°/il,yU2 - 2°/i2) + C3'~(/ii,yU2) eAi2(/il - 2°/ii,yU2 - 2°/i2) 

= -Ai2eAi2(/^i,/^2)- (7.14) 

and since the coefficients Cf^ never vanish or go to infinity on the generic lattice (see 
(I7.12P ). eAi2(/^i ~ 2°/ii,yU2 — 2°yU2) can be obtained for any A12 when the values of exj^^^fii + 
27^1, /i2 + 27^2), eAi2(Aii-2°Aii,/i2 + 27i2), eAi2(yUi + 27ii,/i2-27i2) and ex-,^{fii, H2) are given. 
(See Figure [Ha.) Therefore, if the values of e\^^{^i, 112) are specified on the boundary lines 
/ii = ix\, Hi = fil-2°ni and fi2 = /^s, fJ'2 = /i2-2°jU2 in t he lat tice C±\^^\® C±\^^\, cx-^^if^i, fJ^2) 
can be uniquely determined via the recursion relation (17.141) for all the lattice points with 
/ii < /i* and fi2 < fJ'2- (See Figure [TJb.) 

On the other hand, the LQC operator G/ reduces to the WDW operator in the 
"asymptotic" regime (i.e. l/i/l ^ °fJ^i), and thus the LQC eigenfunction ex{fl) approaches 
the WDW eigenfunction e^(/i) given by (16. 4p for large Therefore, we can use the same 
labels ki to denote the LQC eigenfunctions (namely, Ai,2,3 = ^1,2,2) and similarly decompose 
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e^(/I) as Ckiif^i) ^kjif^^) ^kjyf^s) ■ It follows that Ckjifii) §.kiif^i) fo^' ^ "Z^/ ^^id we have 
the identifications: 



Al2 

and 



167rG'(A;iA;2 + ^), X^^ = 16nG (k^h + ^) , X23 = l6nG [k^k, + (7.15) 



A = uj'\k) = IQttG ( kik2 + kiks + k2k3 + ^ ) • (7.16) 



Like G in the WDW theory, the LQC operator O is not positive definite. For the same 
argument we used in Section IVI Al only those modes with u"^ > are allowed and what 
matters is just the operator 9+, obtained by projecting the action of 6 to the positive 
eigenspace in its spectral decomposition. 

For a given k, we can find e\jj as mentioned above by appropriately specifying the 
boundaries in such a way that e\jj{fii, fij) —>■ ej^jifJ^i) ej^JifJ-j) for \fij\ ^ S> 

However, (16.41) gives Ckiif^i) ^kj'.l-^'i) ~ (4vr)~^|/xi|~"'^^'^'^*'^^ |/i2|~^''^^*'^^ which is not suitable to 
be used for the boundary condition, since one of /ii and fi2 could be small on the boundary 
lines fii = fii = fi* — 2°fii and /i2 = /i2, = Ai2 ~ 27i2 and thus the WDW limit is not 
always valid even if both and fi\ are very big. To reflect the effect of quantum geometry 
for small /ii or fi2, we follow the same strategy a la Thiemann's trick used for the inverse 
triad operators in Section IIVCI and designate the "tamed" function 



1 



X 



7^2 



1 

|/i2 + >2h -|/^2->2| 



(7.17) 



as the value of ex-^^^if^i^ 1^2) on the boundary, which again reduces to CklfJ'i) e^^(/^j) when both 
^ 7^1 and l/iji ^ In practice, we use f l7.17p to specify the boundary condition on the 
lines (/ii = fil, -fi*2 < /^2 < yU2)> (/"i = K-2°yUi, -/i2 < /^2 < IJ-t), (-K < yUi < /i2 = /^2) 
and {—fii < /ii < yUi, /i2 = /^2 ~ 2°yU2) of the lattice for /i^ and /i2 ^ 7^2, then the 

recursion relation fl7.14p uniquely gives the values ex^Jif^i^ (^2) for all the lattice points within 
the range: —fi{ < fii < fi{ and —fi2 < /i2 < fJ'2- The eigenfunctions so obtained are denoted 
as e'l^, ... for the lattices C\^j^\®C\^^\, C\^j^\®C_\^^\ . . . respectively. We are only interested 
in the symmetric sector and have to symmetrize the eigenfunctions to get: 



^'^(/xi, /i2) = ^ (eAit(/^i, ^2) + e^+(/ii, -/i2) + e^it(-/^i' /^a) + e^it(~/^i' "/^s) 
+ ■ ■■ 

+ex~il^i, 1^2) + e^^J/^i, -/^2) + e^J-^i,/i2) + e^^J-//!, -/i2)). (7.18) 

The same procedure is applied for e^^^J^Hi, n^) and e^^^J^fi2, 1^3)- Finally, the LQC eigenfunc- 
tion el'*'' (/i) in the symmetric sector can be determined for a given k by obtaining e^ia , ^Xi2 
and 6^12 separately. 

The eigenfunctions in the symmetric sector of Ti.f^'^ can be normalized such that we have 
the orthonormality relation: 



SM4"(Ar)etf'(,I) *'). (7.19) 
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(where the right hand side is the standard Dirac distribution) and the completeness relation: 

J2 B{jl)ef{jl)^{jl) = ioT\/k, iff ^(^) = (7.20) 

for any symmetric ^I/(/I) G Tif^^. Finally, any symmetric element ^'(/u) of "Hf^^ can be 
expanded as 

^{pi)= f d^k^<{k)ef{pi) (7.21) 



for e'^;\fi)enr\ 

As in the WDW theory, the general symmetric solution to the LQC constraint (17. 5p can 
be written as 

vl/(/7, 0) = / d'k (^+ik) efifl) e-<^ + §„(fc) e^fifl) (7.22) 

where \1/±(A;) are in L^(E:^,d^k). For |^/| > %, (17:22|) reduces to the solution ([63]) to the 
WDW equation. 

If \1/+/_(A;) vanishes, we will say that the solution is negative /positive frequency. Thus, 
every solution to (17.51) admits a natural decomposition into positive- and negative-frequency 
parts, each of which satisfies the first order differential equation in 0: 

-^ = ±^y^^^, (7.23) 

where a/B+ is the square-root of the positive self-adjoint operator 0+ defined via spectral 
decomposition. Therefore, the solution to (17. 5p with the "initial" function \l/±(/i, 0o) = /±(/^) 
is given by: 

^±(/i,0) = e±*v^(*-<^°)^±(/i,0,). (7.24) 



Remark: In the above discussion, we considered a generic e. In the special cases (i.e. e/ = or 
€j = "m), some differences arise because the individual lattices are invariant under the reflection 
/i/ — > —fii (i.e., the lattices and J-'-iejl coincide). For example, let consider the case when 
ei is generic and 62 is special. We have C\^^\ = ^-|e2| = ^^'^ ^Ait ~ ^Ai^" Consequently, the 
symmetrization of (j7.18p should be modified to 

e^^l{fii,H2) = --^(e^+(^i,/X2) + e|+(^i,-//2) + e^+(-/ii,^2) + e0-Aii,-Ai2) 

+eAit(^i'/^2) + e^^^iP-i, -IJ'2) + e^+(-^i,^2) + e;;+(-^i, -^2))- (7.25) 

Furthermore, if both ei and £2 are special, we have e^^^ = e^^^ = e'^^^ = e^^^ and the symmetriza- 
tion is performed as 

e^f^^[^il, ^i2) = --^(e0/ii,^2) + e0/ii,-^2) + 4it(-^i'/^2) +e;[+(-/ii,-^2)). (7.26) 

In the numerical computation, we can focus on the lattice of the special case with e = and 
take advantage of the fact that >C+|£^| = ^-\ej\- Therefore, we restrict ourselves to the positive /^i,2,3 
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and the tedious procedure of symmetrization can be avoided. This way, the numerical computation 
is much more efficient. 

There is a further subtlety when the recursion relation (j7.14p is applied to the lattice L^^ ® L^^ 
with ei = or £2 = 0, because the coefficients may vanish or go to infinity at some lattice 
points and thus more careful consideration is needed. This issue is closely related to the resolution 
of the singularity and the detail of it is studied in Section IVII CI □ 



C. Resolution of the Singularity and Comments on the Planar Collapse 

What is the fate of the classical singularity and planar collapse? In the classical dynamics, 
the Kasner-like solutions (two of positive, the other negative) approach the Kasner-like 
singularity as — oo and the planar collapse as — oo (with > 0), whereas 
the Kasner-unlike solutions (0 < ki,K2,K3 < 1) encounter the Kasner-unlike singularity as 
(f) — — oo and no planar collapse. Since 1 — > 0, according to (13.231) . the singularity (both 
Kasner-like and Kasner-unlike) gives pi,p2,P3 — > 0, which corresponds to the state |0,0,0) 
on Tigi-av thus the point (0,0) on the lattice C^i ® 'C,e2- On the other hand, by (I3.23P 
again, the planar collapse gives (pi,P2,P3) (oo,oo,c)o), which corresponds to the state 
|cxD, oo, oo) on Tigrav ^nd thus the point {oo, oo) on the lattice C^i ® 'C,e2- We first show that 
the state associated with the singularity is decoupled and thus the quantum evolution does 
not break down; later, the behaviors of the planar collapse are commented. 

Unlike the classical theory, the quantum evolution does not stop at the singularity. If one 
begins with a generic lattice, the discrete evolution determined by (17.51) just "jumps" over 
the states with /// = without encountering any subtleties. The subtleties occur only on 
the special lattice as remarked in the end of Section IVII B[ Thus, we revisit the recursion 
relation (I7.14p for the two special lattices: (i) Cei=o ® Ce^^^o and (ii) £ei=o <S> ^^2=0- 

On the special lattice Cei=o ® -^£2^0, the boundary condition given by (17.171) implies that 
6ai2(A''1 = 0,/X2 — 00) = and thus the state \fj,i = 0, /i2 — >■ 00) on TCf^'^Q ^^^^^ (marked o 
in Figure [2]) decouples from the rest part. On the other hand, the values of ex^^^ni, H2) in 
the first and fourth quadrants (marked x in Figure [2la) can be obtained by direct iterating 
(I7.14P while the recursion relation may not work for the values at /xi = —2°^i (points marked 
o in Figure [2lb) since the coefficients 6*^^(0, /i2) blow up by (17.111) . However, this problem 
is merely artificial. To faithfully express the recursion relation, (17. 14p can be rewritten as 

Ct^{fl) exAf^i + 2>i, fi2 + 2%) - C^+{fI) ex,2{f^i - 2°/ii, fi2 + 2%) 
~C+-{fl) ex,2{f^i + 2%,fi2 - 2>2) +C3-"(/i) eAi2(/ii - 2%, fi2 - 2%) 
= -5(/2)Ai2eAi2(/ii,/i2). (7.27) 

Provided 5(0, /i2, /i3)Ai2 = [Note: 5(0,/i2,/^3) = -B(/ii, 0, /is) = -B(/ii, /i2, 0) = by its 
nature], the relation in (17. 7p then gives 

eA,2(2°/il,/i2 + 2>2) = eA,2(2"/il,;U2 -2°/i2) +eAi2(2>l,/i2 + 2°/i2) 

-eA,2(-2>i,/i2 + 2>2), (7.28) 

which can be used to determine the values of eAi2 at the points marked o without diffi- 
culty. Furthermore, if the boundary condition is given symmetric under fii — > — yUi, we then 
have eA,2(27ii,/^2 + 27^2) = ex,2{-2°fii, fi2 + 2712) and (^M) yields eAi2(27^i, /i2 + 27^2) = 
eAi2 (27^1, Ai2 — 27^2), which, as expected, tells that the values at o are the exact mirrors 
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FIG. 2: (a): The lattice vCei=o <^ ^Cej^o- (b): The lattice >Cei=o ® >Ce2=o- The points marked • are 
the boundary lines with given values; in particular, o denotes the boundary points at which the 
value of 6^12 vanishes. The points marked x are the places whose values can be easily obtained, 
while those with o, □ and ® require more careful consideration. (The largeness of \i\ and is 
assumed but not properly presented in this figure.) 

(reflected by the line \i\ = 0) of the lattice points in the line /ii = 2°/ii. Once the values at 
o are all obtained, the iteration can be continued to all unmarked points in Figure Oa. 

Next, let us study the special lattice £^1=0 ® '^^£2=0- Again, we have the states |0, 00) and 
1 00,0) decoupled on the distant boundary (marked o in Figure [Hb). Furthermore, in the 
iterating procedure, the values of 6x^2 can be obtained only to the flrst, second and fourth 
quadrants plus two extra points (all marked x). When the recursion relation is applied for 
the value at the origin (/xi = /i2 = 0, marked (S>), the vanishing of C^~{2°fii, 2712) by (17. lOp 
makes (I7.14p a consistency relation: 

C^^i2%, 2%) eA,2 (4°/ii, A%) - C3-+(2°/ii, 2%) e^,, (0, A%) 
-C3+-(2°/ii,2°/i2) eAi2(4°/ii,0) + Ai2eA,2(2>i,27x2) = (7.29) 

and leaves the value of eAi2(0,0) undetermined. Nevertheless, thanks to fl7.1UI) again, the 
values of eAi2 cit the points marked □ can still be determined without knowing the value at 
®. Moreover, if the Hilbert space Hf^^Q^^^^Q is restricted to the symmetric sector, given the 
symmetric boundary condition, the consistency relation (17.291) is the mirror (reflected by 
the origin) of the relation among the four points marked □ and therefore it is automatically 
satisfled. As a result, the point is completely decoupled and the iteration "jumps" over the 
origin to those points marked □ and consequently to all the unmarked points in Figure [2lb. 

Put all together, it follows that the states associated with the point marked o and in 
Figure [2] are irrelevant. The same arguments can be easily applied for the special lattices 
in (/ii,/i3)- and (/i25 /^3)-planes in the cycUc manner. In the end, combining CaijI/^I! /^2), 
eAi3(/Ui, /is) and eA23(/^2, A'-s) to get e^^\iT), we conclude that the sate |0,0,0) associated 
with the singularity on Tigj-av is completely decoupled from the difference evolution and the 
evolution remains deterministic across the deep Planck regime. Consequently, the classical 
singularity (both Kasner-like and Kasner-unlike) is expected to be resolved in the quantum 
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evolution. (The absence of the Kasner singularity has been shown in [35|] using an effective 
classical Hamiltonian obtained from loop quantization of vacuum Bianchi I model. Here, 
inclusion of the scalar field allows us to prove the absence of the singularity on the exact 
evolution equation as being the emergent time.) 

When the matter sector is included to the model, the directional factor defined as qi := 
P^/IpiI^ for any of the diagonal directions with the constant of motion plays the same 
role as the matter density p := ^p'^/\piP2P3\ does in the isotropic case. The fact that the 
problematic state is decoupled suggests that the directional factor in any diagonal direction 
cannot go to infinity when the quantum discreteness is imposed. Therefore, the classical 
singularity (which gives infinite directional factors) is expected to be resolved and replaced 
by the big bounce in LQC^ However, following the lesson we learned for the isotropic model 
in the critical value of the directional factor can be arbitrarily small by increasing 
since we are now adopting "/^-scheme. Later in Section IVIIIt the dynamics will be refined 
in /i-scheme, by which this problem is supposed to be fixed and the critical value for qj is 
expected to be 0(/i£pf ), independent of Ptp. (This has been shown for the effective dynamics 
in 0.) 

The argument here is essentially the generalization of that used for the isotropic case in |6| . 
The result solidifies the assertion that the key features responsible for the resolution of the 
singularity are robust [2l[ and persist in more complicated cosmological models [sl . However, 
the proof presented here only shows the decoupling of the state in the kinematical Hilbert 
space associated with the classical singularity, but details of the dynamical behaviors of the 
solutions have yet to be investigated. To draw the definitive conclusion that the singularity 
is resolved and replaced by the bounce for those states which are semi-classical at late times, 
we have to figure out the physical Hilbert space, Dirac observables and semi-classical states, 
(which will be the topic of Section IVlIDI) and afterwards a thorough numerical investigation 
can be initiated [111] . 

Remark: The proof in this subsection is rather general, insensitive to the details of the matter 
sector. As far as the matter sector is concerned, the only proviso to make the above argument valid 
is: i?(0, /i2, ^3)Ai2 = -6(^1, 0, /U3)Ai2 = 0, etc. By (|7.15p and ()7.16p . this condition is equivalent 
to: — i?(/2)0(/2) = C'matter(j5) ^ 0, when any of fij goes to zero. That is, C'matter(/^) annihilates 
^(/2, (j)) (0 represents the matter field in general) for /^i = 0, = or /is = 0, which is the case 
for the minimally coupled matter source [41. Il7|. 



Without the discreteness imposed for the inverse triad operator, B{fl) goes to infinity when 
/i/ ^ and the above proviso no longer holds. The finiteness of the eigenvalue of l/^ypj mentioned 
in Section IIV CI therefore plays an important role for the decoupling of the singularity. However, 
the results of both effective dynamics and the detailed numerical analysis in isotropic models tell 
us: it is the "non- locality" of the gravity (i.e, cj replaced by holonomy), not the finiteness of the 
inverse triad operator, that accounts for the occurrence of the big bounce and in fact the bounce 
takes place much earlier before the discreteness correction on the inverse triad operator becomes 
significant. 



^ Note that, unlike the isotropic case, it is the directional factor gi (not the matter density p) that indicates 
the happening of the big bounce. Near the epoch of the singularity, the three directional factors gi may 
approach their critical values at slightly different instants. Therefore, the detailed occurrence of the big 
bounce may happen as many as three times, not just once. 
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Nevertheless, the fact that the problematic state is decoupled kinematically is still important, 
giving the lowest bound at the deep Planck regime, sometime before which the bounce should take 
place. For the semi-classsical states at late times, however, the exact point at which the bounce 
occurs depends on the detailed dynamics and it is very likely to happen much earlier before the 
lowest bound as in the case of the isotropic model. Likewise, in the anisotropic model, we also 
anticipate the same situation, and accordingly expect the directional factor qi as the indication of 



happening of the bounce □ 



Contrary to the resolution of the singularity shown above, the state \oo, 00, 00) on Tigrav 
associated with the planar collapse is not decoupled but persists in the evolution equation. 
This fact is somewhat undesirable since it means one of the length scale factors aj continues 
the vanishing behavior in the forward evolution and thus eventually enters the deep Planck 
length scale without a stop. On the other hand, however, this is expected as the classical 
solutions yield "fiiCj toward the planar collapse and therefore the quantum corrections 
become more and more negligible. (This remains true in u-scheme since /i/C/ — ^ 0. For more 
details of the effectiveness of quantum corrections, see |l2|.) It can also be understood as 
follows. The planar collapse corresponds to the point (cx), 00, 00) in the lattice C^i ®C^2 ®^e3) 
and for the points with large values of pi the lattice discreteness can be ignored. In terms 
of p/, nothing is really collapsing and the classical solution toward the planar collapse is 
perfectly well-behaved if we treat pi as the fundamental variables. 

The result suggests that smallness of pi (not of a/) is the indication for the significance 
of quantum effect and thus also for the occurrence of bounces. While a/ are the length scale 
factors, pi can be interpreted as the area scale factors according to fl3.15p . This seems to 
support the speculation: "area is more fundamental than length in LQG" , although whether 
this is simply a technical issue or reflects some deeper facts is still unclear. (See Section VII. B 



of |36| for some comments on this aspect and [37| for more discussions.) Meanwhile, whereas 
the length operator has been shown to also have a discrete spectrum [s^, the fact that the 
vanishing of the length scale factor is not resolved seems to contradict the discreteness of 
the length spectrum. Whether we have missed some important ingredients when imposing 
the discreteness corrections from the full theory or indeed area is more essential than length 
remains an open question for further investigation. 



D. The Physical Sector and Semi-Classical States 

Results of Section IVII Al and Section IVII Bl show that while the LQC operator O dif- 
fers from the WDW operator G in an interesting way, the structural forms of these two 
Hamiltonian constraint equations are quite the same. Therefore, in the LQC theory, apart 
from the issue of super-selection, we can introduce the Dirac observables and identify the 
physical Hilbert space with the inner product by demanding that the Dirac observables be 
self-adjoint (or by carrying out group averaging method) in the way closely analogous to 
what we did in WDW theory. 

The super-selection sector of the physical Hilbert space '^^phy labeled by e* G [0, 7^1] x 
[0, 7^2] X [0, 7^3] consists of the positive frequency solutions \E'(/i, 0) to (17. 5p with initial data 
\l/(/I, 0o) in the symmetric sector of Hf^^. By f l7.22p . the explicit expression of \l/(/I, 0) in 



33 



terms of the symmetric eigenfunctions ■* reads as 

■q![pi^(l))= ! d^k^{k)ef{fl)e''^'^, (7.30) 

where, as in the WDW theory, uj'^{k) = 167rG(fcifc2 + ^1^3 + ^2^3 + 3/16). The physical inner 
product is given by 

(^l|*2)phy= E Bifl)^,ifI,<f>,)^2{fI,<f>o). (7.31) 

The action of the Dirac observables is independent of e and has the same form as in the 
WDW theory: 



^vl/(;2,0) = e^v^('^-'^°)|/x,|vl>(/2,0„) and p^^^J(fI,^) = -th^^ 



(7.32) 



The kinematical Hilbert space n\°^^^ is non-separable but, because of super-selection, each 

physical sector T^p^y is separable. The eigenvalues of the Dirac observable /i/l,/,^ constitute a 
discrete subset of the real line in each sector. In T^^in^' contrary, the spectrum of pi 

is discrete in a subtler sense: every real value is allowed but the spectrum space has discrete 
topology. Therefore, the more delicate discreteness of the spectrum of pj on H^^^ descends 
to the standard type of discreteness of the Dirac observables on each physical sector Tiphys- 
Furthermore, comparing (17.301) with the WDW semi-classical state given by (I6.17p . by 
analogy, we can construct the LQC semi-classical state at the time (p = 4'o'- 

%„,p;(/^,0o)= / ^ d'k^{k)ef{fl)e^-^''>''"', (7.33) 

where the Fourier amplitude ^{k) is given by 

§(A?) = ^^.(fc)e-^^(^)'^° (7.34) 

with the same ^^(A;) and 0* as defined in f l6.18p and fl6.19p . Since e^'*^(/I) — >• e^(/i) for 
l/x/l 3> the wavefunction given by (17.331) is very close to that in (16.171) for large 
Therefore, the LQC semi-classical state is peaked at the classical trajectory for the region 
far away from the singularity and thus satisfies the pre-classicality [39|, |40|] (at least for one 
asymptotic side of the evolution). 

It has been argued in [41,] that only the zero solution satisfies the pre-classical condition 
for both the vacuum case and the addition of a cosmological constant in the Bianchi I 



model. Inclusion of matter changes the situation as suggested in [41[; here, with being the 
emergent time, the LQC semi-classical sates are constructed and shown to be pre-classical. 

We have shown in Section IVII CI that the quantum evolution is deterministic, but the 
question whether the solution is pre-classical again on the other side across the singularity 
or descends into a quantum foam needs more investigation. Although the thorough numerical 



analysis is awaited to be done [ll[ , the occurrence of the big bounce has been shown in [jj 
and it turns out that the semi-classical states are pre-classical on both asymptotic sides of 
the evolution. (For the issues of pre-classicality and time asymmetry on either side of the 



evolution in the vacuum Bianchi I model, see 42|.) 
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VIII. IMPROVED DYNAMICS (/i-SCHEME) 



Section |V] - Section IVIII complete the analytical investigation for the precursor strategy 
(°/^-scheme). In this section, we repeat the whole analysis for the improved dynamics (/2- 
scheme), which was first proposed and studied for the isotropic model in j^. The prescription 
to endow the fundamental discreteness from the full theory is modified and implemented in 
a more sophisticated way, thus more directly conveying the underlying physics of quantum 
geometry. The corresponding WDW theory is unchanged but recast in Section IVIII Ct 
featuring the use of affine variables (to be defined soon). 



A. Gravitational Constraint Operator 

Fixing the °/i-ambiguity by setting °/i to the finite number can be understood as the 
imprint of the fundamental discreteness of the full theory. In the full theory of LQG, a proper 
regularization procedure is taken to define the Hamiltonian operator: the 3d coordinate space 
is partitioned into small regions of coordinate volume e^. When acting on the diffeomorphism 
equivalence classes, the dependence on e is gone when e is smaller than some particular 
value because of the diffeomorphism invariance. Therefore, the Hamiltonian constraint is 
independent of the regulator in the end, simply because making the partition finer cannot 
change anything below the Planck scale [H, 25. 26l . 



In the symmetry reduced theory, since the diffeomorphism invariance is gauge fixed in the 
beginning, the independence of the regulator is no longer true and we have to fix the regulator 
by hand. However, the strategy used in °/i-scheme is rather heuristic: the holonomy h\^'^ 
in (15. 8p is not a "state" of the Hilbert space and thus setting its "eigenvalue" to the area 
gap is not completely physical. To reflect the fundamental discreteness more rigourously, 
we introduce the improved strategy: "/i-scheme". 

In Cyl of the full theory, the area operator A{S) for a coordinate surface yields quantized 
spectrum when acting on the spin network states |a, j, I). To carry the parallel from the full 
theory to the reduced theory, first, recall that the states Mp. = \ jT) in Cylg are analogous to 
the spin network states J^ajj = in Cyl. Second, note that the area operator A{S) 

in Cyl is replaced by E{Si) '.= E(5/, // = 1) in Cylg defined in f lOD . 

Let 5(^3) be the square surface normal to the z-direction (with respect to °q) with the side 
lengths fl3°Li and fi3°L2 (with respect to °q) in x- and y-directions respectively. According 
to (14.41) . when acting on the state \fT), the area operator associated with S^P-^^ gives 

E(^(^3))|^)=4^y2^^2^3|^). (8.1) 

If S'^'^^) is the smallest area we can shrink to, we should then have Anjip^ fi'^fi^ = ±A or 

l/i^/ial =k:= A/(47r74i) = ^3/2. (8.2) 

In the other two directions, we set S^^'^^ and S^^^^ to the smallest area and have the similar 
results: = 1/^2/^2! = k. Therefore, to faithfully reflect the area gap in the full theory, 

instead of the °^ -scheme used earlier, we should set 



35 



Unlike "/^-scheme, in which "hj are fixed as constant, in the improved strategy, the values fij 
to be set to fix the regulator depend adaptively on the state This is due to the presence 
of the operator pj in (14.41) . 

To figure out what the action of the operator e'^iP-in^^'^ is when acting on we first 

-j-o.. d 

recall that the first line of f l4.17p corresponds to e ^''i /i2, /^s) = ^(/^i =F "fJ'i, f^2, f^s)] 
that is, = in the /T-representation. By the same spirit, we define the affine parameter 

vj such that 

d k d d ^, 



which leads to 



VI = sgn(/i/)^-^ If^if^- (8.5) 
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Changing the coordinates /ii,/i2,Ai3 to fi,f2,f3, we then have 

e^^'^iJj{Hi, 1x2,1^3) = e^^i!{vi,V2,V3) = ij{vi =f l,V2,V3), (8.6) 

or equivalently 

e±5/^^"Vl,Ai2,/^3) = \Vi±l,V2,V3). (8.7) 

In terms of the new variables, that is 





It/) 











|t;i,t;2,t^3±l). (8.8) 

In /i-scheme, the new operators corresponding to fl5.9p and flS.lOp are obtained by re- 
placing 7^/ with fij. Consequently, when acting on the states labeled by 1^1,^2,^3, (15.110 is 
modified to 

Cjil\vi,V2,V3) 



167r73£2j^^^2/i3 



x<! (IK 1+2, f2+2, 113+1 l+2,t'2+2,l)3- 



-ll + IKl ,f2, 1^3 + 1 Kl,i;2,«3 — 1 

|K)i-2,l'2+2,l)3 + l ~ K)i-2,t)2+2,t>3-l ,f2,f3 + l ^1,'1'2,1'3 — 1 

\]\vi- 2,V2 + 2,t;3) 

~ (^l^fi+2,i;2-2,t)3+l ^ K)i+2,t)2-2,t>3-l| + |K>i,t)2,D3+l ~ ^vx,V2,V3-l\^ \Vl + 2, f2 — 2, fs) 
+ ^|K;i-2,t>2-2,/i3+l — K)i-2,i>2-2,i3-l| + |K'i,i)2,i'3+l ^ K)i,i;2,f3-l|) 1^1 ^ 2,f2 — 2,^3) 

+ cyclic terms > (8.9) 



In the ordinary quantum mechanics, P — and e''°'^^^'ip{x) — ijj{x + a) but we cannot simply take 

^if{x)p/h _ g/(2:)^ because it is not unitary in generaL Here, on the contrary, since Cyl^ is in the "polymer 
representation" and the right hand side of (|4.7p is the Kronecker delta (not the Dirac distribution), the 
same problem does not occur. In fact, 'H^{^ = L^(KgQjjj., d^^Bohr) = ^^(l^Bohri '^'^"Bohr) and e^^^(^''^*^' is 
unitary. 
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and fl5.12p is modified to 



167r7¥|iyUi/i2/i3 

,1)2+2, 113-1 I 1^1 + 2,^2 + 2,^3) 

+ |Ki,-!;2+2,i)3+l — K;i,-U2+2,-!;3-l|^ ^1 " 2, t>2 + 2, fs) 

K, I |i;i + 2,t;2-2,i;3) 

~l~ I Kll, 1)2 — 2,1)3 + 1 K)i,i;2 — 2,1)3 — 1 



-(IK 
-(IK 

+ (|K 

+ cyclic terms|, 
where the eigenvalues of the volume operator are given by 



'111—2,112,1)3 + 1 Kll — 2,D2,1'3 — 1 I 
'1)1+2,112,1)3 + 1 K)i+2,l)2,f3 — 1 I 
'di— 2,112,1)3 + 1 K)i — 2,1)2,1)3-1 I 



V,, 



Dl,D2,l'3 



)|t;i-2,t;2-2,t;3) 

(8.10) 

(8.11) 



In terms of the new notation, it is equivalent to 



C'Svl^l,^^2,^^3) 



vr^pi 



272^l/i2/U3 



where 



3Vk 



C++iv) \v, + 2,V2 + 2, vs) - C,+iv) \v, -2,V2 + 2, V3) 

-Ctiv) \v, + 2,V2- 2,vs) + C^^iv) -2,V2- 2,V3) 
+ cychc terms|, (8-12) 



\{v, + 2)v2iv3 + 1)1^/3 -\iv, + 2)v2ivs - 1)1^/3 



3Vk / 



+ |^^i(^^2 + 2){v, + 1)|V3 _ |^^(^, + 2){v, - 1)1^/3 
\{v, - 2)v2{v, + 1)1^/3 - \{v, - 2)v2{v, - 1)1^/3 



3v^ 



+ |t;i(t;2 - 2){v, + 1)1^/=^ - \v,{v2 - 2){v, - 
\{v, + 2)v2{v, + 1)1^/=^ - \{v, + 2)t;2(t;3 - 



+ 



|^^i(^^2 - 2)(t;3 + 1)1^/^ - W{V2 - 2){vs - 



\{v^ - 2)V2{V3 + 1)1'/=^ -\{v,- 2)V2{V3 - 1)| 



1/3 



+ 



\vi{v2 + 2)(t;3 + 1)1^/3 - \v,{v2 + 2)(t;3 - 1)|^/^|) 



(8.13) 



and similar for C^^ and Cf * in the cyclic manner. 



By inspection, it is clear that CgrL and Cgrav commute with the parity operators 1112,3: 



[c^^l, Ui] = 0, [^(f),, n,] = 0, 
37 



(8.14) 



where Ili\v) = |— "Wi, f2, "^3) (or equivalently, Ilil/i) = |— /ii, /i2, /is)) and so on, which flip the 
orientation of the triad. 

Note that, unhke CgJ^l and C*grav, because of the presence of nonconstant j2i in fl8.9p 
and (IS.lOp . the operators CgiL and Cgfiv are not self-adjoint in Tigrav; unless further re- 
ordering is performed. Nevertheless, the operator to be defined in f l8.2ip is self-adjoint in 
L'^iM.^^^j., B{v) d^VBohr) with proper measure B{v) chosen. Therefore, the procedures used 
in °/i-scheme to construct the physical sector can still be applied for the improved scheme. 
(Also see the remark in Appendix IB 21 for the issue of re-ordering in the isotropic case.) 

Remark: To impose the "adaptive" regularization, a seemingly direct but in fact misconceived 
method would be to set S^^^^'^^^ (the rectangular surface normal to the z-direction with the side 
lengths p-3°Li and fL3°L2 in x- and y-directions) to the smallest area A and so on for the other two 

directions. Instead of ([83]), this method gives fli = \J}^\, ^2 = \J jj^^, and //s = \[}^\, which 
mingle the three diagonal directions and complicate the Hamiltonian constraint (in particular, the 
evolution constraint is no longer a difference equation). A more careful consideration urges that 
we should use the "square" surface S^'^^^ rather than the generic "rectangular" surface S^^'^'^^^ to 
reflect the discreteness imprinted from the full theory. In the full theory, the discretized area is 
associated with the edge of the spin network which penetrates the surface and carries the su(2) 
label j. Therefore, in the symmetry reduced theory, to properly impose the discreteness, the role 
of j is replaced by the single label fij in the normal direction to describe the area of the associated 
surface.^^ □ 



B. Inverse Triad Operators and the Total Hamiltonian Constraint 

In /2-scheme, by the same spirit as discussed in the previous subsection, the inverse 
triad operators defined in (14.151) should be changed as well by replacing °hj with fli, and 
correspondingly, when acting on the states, fl4.16p is modified as 



\v) 



where Lj is the eigenvalue of the length operators: 



\vi + 11 



\vj - ih 



(4717) 



1/2 



1/3 



(8.15) 



(8.16) 



Note that as in "/^-scheme, the inverse triad operators in /i-scheme commute with pj and are 
self-adjoint in TC^^-^ and bounded above as well. The upper bound is given at vj = ±1: 



bli = 2^/'(4vr7A:)-^/^^p/. 



(8.17) 



The author thanks Ashtekar and Pawlowski for clarifying this confusion. 
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Furthermore, the eigenvalue of 1/ ^^/pj goes to zero when fij ^ and this fact is important 
for the resolution of the singularity as discussed in Section IIV CI 

Since = SirGp'^/ \/\piP2P3\ in fl3.22p . the corresponding operator acting on the 



states of mt''^ 



''kin 



nL.. ® yields 



''grav 



Cf\v) 



(4717)2 £3 



■B{v)\v) (»87tGpI\ 



(8.18) 



PI 



where 



B{v) 



3Vk 



2 /il/X2/^3 



\vi + 11^/^ -\vi-l 



1/3 





X 


1^^3 + 1 


1/3 _ 


\V2 


- 1 


1/3 


9 


1 |l/3 


l^^l + 1 


1/3 _ 


\vi 


- 1 


1/3 




X 


1^^3 + 1 


1/3 _ 


\V2 


- 1 


1/3 



\V2 + 1 



V2 



1/3 



\V2 - 1 



1/3 



1 1/3 



\V2 - 1 



1/3 



The total Hamiltonian constraint Cg^v + C'^f^ = then gives 

87rG B{v) Pi \v,^) + (47r7)i4iC'Sv ^ 



0. 



(8.19) 



(8.20) 



Again, we choose the standard Schrodinger representation for but the "polymer repre- 



sentation" for V (or equivalently for /I); that is Ti. 
this representation, p, 



total 
kin 



Bohr) '^^'^^Bohr) < 

and the total Hamiltonian constraint read as 



dcf)) . In 



-C+-{v) ^{v, + 2,V2- 2, vs, 0) + C^-{v) ^{vi -2,V2- 2, v^, 
+ cyclic terms I 

{Q3{VI,V2) + Q2{VI,V3) + Qi{v2,V3))^{v,^) =: -e{v) ^ {v, (f)) , 



(8.21) 



where we define 



B{v) := fiifi2^i3B{v) 



\vi + l 



1/3 



\vi-l 



1/3 



1^^2 + 1 



1/3 



1^2-1 



1 1/3 



1^^3 + 1 



1/3 



\V2-I 



1 1/3 



(8.22) 



By inspection, we can easily see that G and each of 61,2,3 individually are self-adjoint on 
L\Ri^^^,B{v)d^VBoi.r). 

Note : The total Hamiltonian constraint in /i-scheme for the isotropic case can be obtained by 
formally "translating" the anisotropic results found in this section. This is shown in Appendix IB 2\ 
where the different factor orderings for C^^^ and Cgriv used in this paper and in 1^] are also 
remarked. □ 
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C. WDW Theory (Revisited) 

When the improved dynamics (/i-scheme) is adopted in LQC, the corresponding WDW 
theory remains the same. In order to compare the WDW theory with LQC in the new 
scheme, we recast the results in Section I VI I in terms of the affine variables vi defined in 
fl8.5p . This change of variables is straightforward and we only list the results. 

The total WDW Hamiltonian constraint in fl6.3p becomes 



_ i9 /, 1, 2\ d , ,1 d /. 1. 2\ d 

X I l^^2pT^ K^2r K^i r ^ — \- \vi\^^r~ Fir ^"2^ tt- 

OV2 V / ovi ovi V / 0V2 



-cyclic terms 



o ^rn/-NT-i2 f SVky rsgn(fif2) f d . .id d . .id 
= 8.G|B(„)| '-(^— j {_i^(^_K„,|._ + _|,,„,|._ 

+cyclic terms I 

=: -iasi^i,V2) + Q2ivi,V3) + Q^iv^,V3))^ =: ~Qiv)^, (8.23) 
where B_{v) = B_{fl) = (3v^/2)^^|fif2f3p^^^ and we also define 

B(v) := (? J (^"j \:,v.v..r"m = (5 J (^) l«.«2»3|-'^ (8.24) 



which gives 

5(/i) d^fi = B{v) d^v. (8.25) 

Note that in LQC theory, B{v) and B{v\ defined in fl8.19p and fl8.22p give the expected 
asymptotic behaviors: B{v) B_{v) and B{v) B_{v) when |fi,2,3| ^ 1- 

The operators 9 and each of 6^^ 2 3 individually are self-adjoint on L^(]R^, B_{v) d^v). The 
(symmetric) eigenfunctions of 6 in (16. 4p now reads as 

f3v^/2r^^^^'''^'^'^ 1 

e^{v) = e^(/2) = ^ ' ^^^^^^ \VlV2V3\-^ ^ipilnM giffeln^^l gi|fc3lnb3|_ (8_26) 



The corresponding eigenvalues are unchanged and given by (16. 50 . The orthonormality rela- 
tion (16. 6p is rewritten as 



dhB{v)e,;{v)e^,{v) = 6%k, k'), (8.27) 
and the completeness relation (16. 7p gives 

j d\B{v) e^{v) ^(tT) = for V^, iff ^(tJ) = (8.28) 
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for any ^{v) G L|(R3, d^v). 

The complete set of Dirac observables is comprised of and v\^^ {v at some particular 
instant (j) = (j)o). The representation of the Dirac observables on TC^^ is given by (cf. 
(EIE])): 

=e'^^'^"'^°V/l^(^,0o) and p^<iliv, (j)) = hy^^{v, (I)). (8.29) 

Both ptf) and vl^p^ are self- adjoint on the physical Hilbert space 'H^^ endowed with the 
inner product (cf. ( 16.15^ ): 

(vI/i|M/2)phy = jdhB{v)^,{v,<p,)<il2iv,<Po). (8.30) 

The semi-classical (coherent) states are the same as discussed in Section IVI CI The 
expectation value of {71 is given by (cf. (16.211) ): 



2 



3/2, , . 2 



{viUo) = sgn(/i/)^^( ) ^sgn{fii)^^\{fii\^^)f'^ 

^ sgn{f,j)^\f,%f (8.31) 

and the dispersion square of vl^p^ is (cf. f l6.22p ): 

The approximations used here are valid if the semi-classical states are sharply peaked at the 
classical values. 

The WDW theory recast in terms of the affine variables gives us the asymptotic behavior 
(|'?^i,2,3| ^ 1) of the general solutions and the physical sector of LQC, which will be discussed 
in the next subsection. 



D. Emergent Time, General Solutions and the Physical Sector 

Just like the Hamiltonian equation (17. 5p in °/i-scheme, the improved Hamiltonian equation 
(I8.2ip in /i-scheme has the similar structure. All the results studied in Section IVllBI - 
Section rVlIDI can be easily reproduced and modified for the /i-scheme in the obvious manner 
(using the affine variables vj). Again, fl8.2ip can be regarded as an evolution equation with 
respect to the emergent time and the 0- independent operator Q{v) is a difference operator. 

The space of physical states is naturally divided into sectors, each of which is pre- 
served by the evolution and by the action of the Dirac observables. This again implies 
the super-selection for different sectors. Let C\^j\ (|e/| < 1) be the "lattice" of points 
{|e/| -|- 2n;n G Z} on the f/-axis, be the "lattice" of points { — le/l + 2n;n G Z} 

and Cej = jC-iej] [jC-\ej\- Also let 'H^^^^^ _|_|^^| ^^^^^ (8 of them) and TCf^^ denote the physical 

subspaces of L^(R|qJjj., i?('i7) (i'^t>Bohr) with states whose support is restricted to the lattices 
®C±\e2\ ®'(^±\t3\ and Cei ®Ce2 ® respectively. Each of these subspaces is mapped to 
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itself by O. Furthermore, since CgrL is self-adjoint on H^^^ = L'^{M!^^^^,d^VBohi:), it follows 
that is self-adjoint on all these physical Hilbert subspaces. 

As in the 7i-scheme, since C\^j\ and are mapped to each other by the operator IT/, 

we restrict ourselves to the symmetric sector (i.e. the subspace of Tif'^'' in which ^(v) = 

For a generic e (i.e., none of e/ equals to or 1), we can solve for the eigenvalue equation 
6 e\{v) = (01 + 02 + 03) ex{v) = A e\{v) on each of the 8 Hilbert space 'H^^^^^ _|_|^^| _|_|^^|. Let 
ex{vi,V2,v-i) = ex^{vi) ex2{v2) ex^^vs) and solve (cf. (|7.13p): 

03eAi2(^^i,^^2) = Ai2eAi2(t;i,t;2), ioT ex^2ivi,V2) = ex^ivi)ex2iv2), 
02eAi3(fi,f3) = Ai3eAi3(fi,f3), for ex^.,ivi,V3) = ex^ivi)ex:,ivs), 
0ieA23(^^2,^^3) = A23eA23(^^2,^^3), for eA23 (t^2, ^^3) = eA2 (^^2) CAg (t^s) , (8-33) 

with A = A12 + A13 + A23. (Note that 0/ is independent of vj.) Each of fl8.33p can be solved 
recursively for any arbitrary Ajj provided that the appropriate boundary condition is given. 
For example, 03 eAi2('?^i5 '^^2) = A12 eAi2(^i5 ^^2) reads as (cf. (17.141) ): 

C^3^+(^^i, ^^2) ex,,ivi + 2,V2 + 2)- C^+{vi,V2) ex,,ivi -2,V2 + 2) 
-C^'{vi,V2) 6x^2(^1 + 2,f2 - 22) + C;^~{vi,V2) 6x^2(^1 -2,V2- 2) 

= -Ai2eAi2(t^i,t^2), (8.34) 
where we define 

Cf^{vj,VK) := ^[B{v)]-'Cf^{v), (8.35) 

which is independent of /i/ and has the same properties as Cf^{^j, ^k) does in f lT.lOp - 
f l7.12p . Since the coefficients never vanish or go to infinity on the generic lattice, 

eAi2('^i — 2,112 — 2) can be obtained for any A12 when the values of exi2{'Vi + 2,^2 + 2), 
^Ai2 ("^1 ~ 2, t>2 + 2), eAi2 {vi + 2, f 2 — 2) and eAi2 ("^^i, '^^2) are given. As a result, if the values of 
eAi2('^^i; "^2) are specified on the boundary lines Vi = v^, Vi = v* — 2 and V2 = ^2, f 2 = — 2 
in the lattice C±\ej^\^C±\e2\, exi2{vi,V2) can be uniquely determined via the recursion relation 
fl8.34p for all the lattice points with vi < and V2 < 

On the other hand, the LQC operator 0/ reduces to the WDW operator 0j in the asymp- 
totic regime (i.e. \vi\ ^ 1), and thus the LQC eigenfunction ex{v) approaches the WDW 
eigenfunction e^('y) given by fl8.26p for large Therefore, we can use the same labels 

ki to denote the LQC eigenfunctions (i.e., Ai^2,3 = ^1,2,2) and similarly decompose e^('y) as 
Cfci('^^i) 6/02(^2) e^3(^^3)- It follows that Ckjivi) — >• e^ii'^i) ^ 1 and we have the identi- 

fications (IT.lSp and fl7.16p unchanged. Again, the LQC operator is not positive definite 
and only those modes with uj"^ > are allowed. What matters is the operator 0_|_, obtained 
by projecting the action of to the positive eigenspace in its spectral decomposition. 

For a given k, exjj can be obtained as mentioned above by appropriately specifying 
the boundaries in such a way that exjj{vi,vj) — > ei.^{vj) Cj^j^vj) for |t>/| ^ 1, \vj\ ^ 1. 
To be used for the boundary condition, however, the asymptotic function Ckii'^i) ^hi'^^'i) = 

(47r)-2(3|l)-|+f (fci+fc2)|^^|~|+f'=i|y2r^^^''' given by fl8:26|) should be further "tamed". 
Following the same spirit used for °/i-scheme, the value of eAi2(/^i) /^2) on the boundary 
is given by (cf. fl7.17p ): 



(3v^/2)"^+-^'^^^"^ /3 , ^,1 , .,1 V-'^Y3 



(47r)2 ^2 



X 



|t;2 + lh -|f2-l| 



(8.36) 
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which again reduces to Ckji'^j) ^kii'^j) when both \vi\ ^ 1 and \vj\ ^ 1. In practice, 
we use ( I8.36P to specify the boundary condition on the hues (fi = v^, —V2 < f 2 < '^^2)5 
{vi = f ^ — 2, —V2 < V2 < V2), (— fi < f 1 < f^, V2 = V2) and (— < f 1 < v^, V2 = V2 — 2) of 
the lattice for 1 and f 2 ^ 1, then the recursion relation fl8.34l) uniquely gives the values 
ex-^2ivi,V2) for all the lattice points within the range —v^ < vi < v* and —V2 < V2 < V2- 

Similarly, all the other results and remarks in Section IVII B| - Section IVII Dl can be easily 
modified for the /i-scheme by formally taking the replacements: ni — >■ vj and B{fl) — >■ B{v). 
In particular, the orthonormality relation for the eigenfunctions in the symmetric sector of 
nf" is given by (cf. fl7J9|) ): 

J2 m ef{v) ef^v) = 6\k, k') (8.37) 

and the completeness relation reads as (cf. (17.201) ): 

B{v) ef{v)-^{v) = for VA?, iff ^(tJ) = (8.38) 

for any symmetric "^{v) G Ti?^^ . The general symmetric solution to the improved LQC 
constraint (18.211) can be written as (cf. (17.221) ): 



d^k (^+{k) ei'\v) e'^'^ + ^^{k) ef{v) e-'^'^\ (8.39) 



t^2(A,.)>o 



where ^±{k) are in L^(M^, d^/c). The solution to fl8.2ip with the "initial" function 
'^±{vAo) = f±{v)i^ (cf. (EMD): 



m^{v,(t>) = e±Ve+(*-<^°)^i(^,0„), (8.40) 

where a/0+ is the square-root of the positive self-adjoint operator 6+ defined via spectral 
decomposition. 

Finally, the super-selection sector of the physical Hilbert space Hp^y labeled by e e 
[0, 1] X [0, 1] X [0, 1] consists of the positive frequency solutions \1/('U, 0) to (I8.2ip with initial 
data '^{v,(l)o) in the symmetric sector of Tif^^. By (18.391) . the exphcit expression in terms 

of the symmetric eigenfunctions e'^^ is given by (cf. (17.301) ): 

*(v,0)=/ d^k^!{k)ef{v)e''^'t', (8.41) 

where uj'^(k) = 16rrG{kik2 + kiks + /C2/C3 -|- 3/16). The physical inner product on the physical 
sector Hphys is given by (cf. (17.311) ): 

(*i|^2)phy= Yl B{v)^i{v,<Po)'^2{v,<Po). (8.42) 

The action of the Dirac observables is independent of e and has the same form as in the 
WDW theory (cf. (TTm ): 

^^^(t;,0) = e*v^('^-<^°)|t;/|^(t;,0o) and ^^^({7,0) = -?^^^^^. (8.43) 
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Just like the case in "/^-scheme, because of super-selection, each physical sector Hp^y is 
separable while the kinematical Hilbert space is not. The standard type of discreteness 

of the Dirac observables f/|<^o on each physical sector is descended from the more delicate 
discreteness of the spectrum of vi on 71^^ = L'^(M.^^^^, (i^fBohr)- 

Finally, it is easy to check that the coefficients Cf^{v) and Cf^{vj,VK) have the prop- 
erties exactly the same as those of Cf^{p) and Cf^{fij^ fix) listed in (17.71) . (17. 8p and 
(I7A0|) - (I7:T2|) . Together with the fa ct that B{0,V2,V3) = B{vi,0,V3) = B{vi,V2,0) = 0, 
the argument used in Section IVII CI can be applied to prove the resolution of the classical 
singularity in /i-scheme and to suggest the occurrence of the big bounce. Moreover, by the 
lesson learned in [8|, in /i-scheme, the bounces are expected to happen when the directional 
factor Qj := p^/lpif enters the Planck regime of 0{hipf), independent of the values of p^. 
Meanwhile, the planar collapse is not resolved but persists in /i-scheme, yielding the same 
comments as made in the end of Section fVII C[ 



IX. DISCUSSION 

The analytical investigation in this paper is the direct extension of jl| and ^ , developing a 
comprehensive LQC framework for the anisotropic Bianchi I model. Both the precursor 
) and the improved (ft-) strategies are successfully deployed. The key features responsible 
for the desired results remain the same: (i) the kinematical framework of LQC forces us 
to replace the connection variables by the holonomies, and (ii) the underlying quantum 
geometry gives rise to the discreteness in the spectrum of the triad operators, as the imprint 
of the fundamental "area gap" in the full theory. Furthermore, following the /i-scheme newly 
proposed in [8|, a more sophisticated prescription is introduced to impose the discreteness, 
thereby capturing the underlying physics of quantum geometry in a more sensible way. 

The success of this extension affirms the robustness of the methods in [sl and thus 
enlarges their domain of validity. It also shows the fact that the singularity persists in the 
WDW theory is rather generic, even in the less symmetric model. The main difference 
with anisotropy is that the operator B is no longer positive definite. Nevertheless, the de- 
tailed analysis shows that the only relevant part is 6+ (the projection of 9 on its positive 
eigenspace) and therefore the strategies can be readily carried through without a big modi- 
fication. Consequently, the expected results for the Bianchi I model are established: (i) the 
scalar field again serves as an internal clock and is treated as emergent time; (ii) the total 
Hamiltonian constraint is derived in both °/i- and /i-schemes and gives the evolution as a 
difference equation; and (iii) the physical Hilbert space, Dirac observables and semi-classical 
states are constructed rigorously. 

On the other hand, the variety of anisotropy does give new features, notably distinct from 
those in isotropic models. Unlike the isotropic case, the Bianchi I model with a massless 
scalar field admits both Kasner-like and Kasner-unlike solutions in the classical dynamics. In 
the quantum theory, we anticipate that the singularity (both Kasner-like and Kasner-unlike) 
is resolved. This is evidenced by the observation: in the kinematical Hilbert space, the state 
associated with the classical singularity is completely decoupled from the other states in the 
difference evolution equation, and therefore the evolution is deterministic across the deep 
Planck regime. 

Meanwhile, we have the somewhat surprising observation that the planar collapse for 
Kasner-like solutions is not resolved and thus the vanishing behavior of the length scale 
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factors aj is not prevented by the quantum effect. This suggests that either area is more 
fundamental than length in LQG or some important ingredient is missing when we impose 
the discreteness correction from the full theory. Either case requires further studies. 

To actually figure out the behaviors of the quantum evolution, the analytical investigation 
alone is not e nough. A complete numerical analysis or at least the study of effective dynamics 
is in demand [ll|, ll2[]. We expect that the big bounce occurs whenever one of the area scales 



Pj undergoes the vanishing behavior. It may happens as many as three times. Furthermore, 
in the improved dynamics (/2-scheme), the bounce takes place at the moment when the 
directional factor gi in any of the three diagonal directions comes close to the critical value 
in the Planck regime, regardless of the values of p<^. 

Other possible extensions have been discussed and commented in Discussion of In 
particular, in order to understand the BKL scenario in the context of LQG, the next step 
would be to explore Bianchi IX models based on the same formulation. Meanwhile, as 
both the isotropic and Bianchi I models have been analyzed in great detail, another feasible 
research direction is to study the emergence of isotropy in the arena of Bianchi I models, 
which could shed some light on the symmetry reduction for more complicated situations. 
As remarked in the end of Appendix IB 1\ this issue is not trivial and demands more detailed 
investigation. 
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APPENDIX A: WDW SEMI-CLASSICAL (COHERENT) STATES 

In this appendix, we study more details about the WDW semi-classical states described 
in Section IVI GI We first look closer at the semi-classical state to better understand its 
structure and then compute the physical quantities corresponding to the Dirac observables. 

The semi-classical states in WDW theory are constructed in fl6.17p with ^^^^(fc) sharply 
peaked at = /Ck. Since \E'k(/c) is sharply distributed, only those k close to JCk contribute. 
For /Ck0 ^ 1, u!'^{k) ^ 16irG{kik2 + kik^ + ^2^3) ~ SttGK.'^k^ in the vicinity of k = Kk. 

As a result, the physical condition cj^ > is automatically satisfied for the relevant k. 
Therefore, when a sharply-peaked semi-classical state is concerned, the restriction cj^ > 
for the integral range in (16.171) can be ignored. We will make use of this fact to simplify the 
following calculations. 

With the inner product defined in fl6.15p . we have 



dkidk2dk3 e ^"'t e e ^"i =1, (Al) 



(27r)3/Via2a3 



00 fCO POO (fci-KKi)2 (fe2-KK2)^ (fc3-KK3)2 



—00^ —00 
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where the identity (I6.6P is used and the restriction > is released. The overall factors in 
fl6.18p has been chosen to make the physical norm of the state normalized. 

Before computing the expectation value and dispersion for fl\ij>^, first note that by 
e°-d^ f{x) = f{x + a) and (16. 4p we have 

e"^"^e,,(/7) = \^^JreM (A2) 

and consequently 

J d%B{f^) Ufl) \f^ire^,{fl) = e-^'^S'ik, k'). (A3) 
The derivatives of 6^{k, k') can be handled with integral by part. This then leads to 



d^k Sk' 6%k, k') ^!^{k') e-^'^('=')(0o-0*) g*M7 (^!f^{k) e^-W('^°-<^S) j . (A4) 

Applying e""^ f{x) = f{x + a) again, we have (taking / = 1) 

e'A (^^^{k) e*^W(-^°-<^S) j = ^^{ki + i, k2, k^) ei^{k^+iMM){^o-ro) 



e (2<Tij^ g 



= e 

This follows 



d^ke e VV'»i*2+fei*3+*2'=3/'^° ^.(A;)^^(A;) 

~ e(2^e ^ ^ d^k^ft{k)^i^{k) 



1^ V^(^i^)(</'*-0o) 



/^i<^„, (A6) 



where in the second line, we simply assign the peaked value to the integrand. This is the 
result in dOT]) . 

Similarly, taking a = 2 in (]A2p and following the same procedure, we have 



4 2^8^ 



= e(^e ^ '^^ ^'"^^ =e^(/it^j2, (A7) 

which leads to the variance in (16.221) . 

Next, to compute the expectation value for p^, we recall the action of p</, defined in (I6.14p 
and use the fact a/0_|_ eg(/^) = ^{k) ej:{jT). It follows 

{p^) = f d^khuj{k)^^{k)^R{k)^huj{lCK) f d^k^S)^S) 

Juj^{k)>0 J 

= hVSnGJCK^ = pI, (A8) 
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where again we ignore the restriction for the integral range and directly put the peaked value 
to the integrand. This is the result in fl6.23p . 
Following the same procedure, we have 



{pD = [ d^kh^u:\k)^S)^S) 



'<^2(fc)>0 

d k{kik2 + kik^ + k2k^) e 



(27r)3/Via2a3. 

= lQ'KGh^lC{Kin2 + Kl/t3 + /t2/t3) = Pf, (A9) 

which leads to A'^{p^) = — (p</,)^ = 0. However, the vanishing of A{p^) is artificial, 
resulting merely from the approximations we have used in ( lASp and ( ]A9p . which neglect the 
variation in A(p<^) in the order of order to recover the correct estimate of 

A(p0) for the desired order while keep the integration tractable, we assume {p^) ~ (p^)^''^ 

(based on the assumption that is sharply peaked) and compute (pi) instead: 



(pj) = / d'kh'u;\k)^,{k)^,{k) 

- / d^k{kik2 + kik3 + k2k3f e 



(27r)3/Via2a3, 

(167rG)2;i^|/C2(a2(/t2 + + al{K^ + K:,f + al{K, + K2f) 

+ {'^1^2 + '^1^3 + ^2^3) + K,'^{f^lf^2 + f^lf^3 + '«2'«3)^| 



4 



icTy2+alal + ay,)y (AlO) 



This then gives 

-^((^ni (^") 

and therefore 

A'(P0) = {pD - {P<t>? ^ T^{^!i^2 + + ali^i + + a|(«:i + /^2)'), (A12) 
which is the result shown in (16.241) . 



APPENDIX B: REMARKS ON THE ISOTROPIC MODEL 
1. Precursor Strategy (/io-Scheme) 

In this subsection, we discuss how to read off the gravitational constraint operator for the 
isotropic case from the results obtained in Section IV Al and Section IVBI for the anisotropic 
cases. In the end, we comment on the symmetry reduction from anisotropy to isotropy. 
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In the isotropic model, the classical phase space is given by the variables c and p with 
the fundamental Poisson bracket (note the difference of the factor 3): 



{c,p} 



(Bl) 



In quantum kinematics, the eigenstates of p are the basis states 



PI 



(B2) 



The basis states are also the eigenstates of the volume operator V: 

3 

V\fi) = V,\fi) := (^\^,\Y ei^\^). 



(B3) 



To reproduce the isotropic result from the anisotropic Bianchi I models, we can simply 
take the replacement: 

c^c^c^ — > c, pi,P2,P3 — 'p- (B4) 

In LQC, correspondingly, we take 



l/ii + a°fii,fi2 + (3°ii2,rn^ + — ^ |/i + (a + /3 + 7)/io) 



and 



^l+«°/xi,/i2+/9°At2,A»3+7°M3 ^+(a+/9+7)/io5 

where in the end we also impose 

= 3°fii = 3k^ 3V3/2. 



(B5) 
(B6) 

(B7) 



The replacement rules then "translate" ([52]) O, flSTTOD and ([EZD flSAT]) fl512D to their 
counterparts in the isotropic model: 



_ 24i sgn(p) 



sin^(/ioc) ( sin(^^)V^cos(^^) — cos(^^)V sin(^^ 



/^(Mo)' _ 1 [/^(Ato) _|_ /^(Ato)t 
^grav 2 grav,o ~ ^grav,o 



A(Mo) _ 24i sgn(p) 



(B8) 



(B9) 



sin(M (sin(^)V^cos(^) - cos(^)V sin(^) ) sin(M (BIO) 



and 



(7^''°) \a) 



^77734^/^3 



K+f.. - V/^J (lAi + 4/io) - 2|/i) + l/i - 4/io)) 



(Bll) 



^grav I A'') 



1677734^/13 

+ (| - Vs^J + I W - Wl) I/" - - 4|V^+^„ - Wll/i)}, (B12) 
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6 r 

-(I - Wl + I W - VsmJ)!-")}- (B13) 

The last one with /Xq = 3/c is the master equation used in 



Remark: According to (j5.12p . the subspace spanned by /i2, /U3) with /Ui = /i2 = is not 
invariant under Cgrav- Hence, by restricting the initial states to the "isotropic sector", the total 
Hamiltonian constraint Cgrav + =0 does not yield the solution which lays on the isotropic 
sector for all the time of the evolution. To systematically study how the isotropy appears in the 
anisotropic framework, more conceptual considerations and technical tools for a better notion of 



symmetry reduction such as developed in 43, l4J] may be needed for the future research. □ 



2. Improved Strategy (//-Scheme) 

In this subsection, we follow the ideas used in Section IVIII Al and Section IVIII Bl to 
construct the total Hamiltonian constraint in /i-scheme for the isotropic case. In the end, 
the issue of the different factor orderings for and Cgfiv used in this paper and in 

are also remarked. 

In the isotropic model, the area operator in (14.41) is replaced by 

E{SJ)=pVo~~'Asj, (B14) 

with Asj the area of S measured by °qab times the factor /. Consider S^^^ to be a square 
surface with side length fiVo (measured by °qab), which gives 

E{s^-^^)\^,)=p^'\^,) = (B15) 

When the /i-scheme is applied to the isotropic case, again, we require the smallest area to 
be 47r7£pi/i^yu/3 = ±A and thus we set 



/X = vWH- (B16) 
Following the same procedure for the anisotropic case, we solve the affine parameter 



and it leads to 

e 

or equivalently 

Therefore, (1B13P is altered to 



= e^9^ip{v) =ip{vTl) (B18) 
e^-^f^'lfi) = \v±l). (B19) 



Cit \v) = C^v)\v + 4) + C~{v)\v - 4) - C%v)\v) [, (B20) 



r/3 



72/1' 
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where 



C-{v) 
C\v) :-- 



3^ 


/3A; 




2 


3^ 






2 


3^ 


/3fc 



v + ?>\ -\v + l\ , 
v-l\-\v-2>\ , 
v + ?,\-\v + 1|| + \\v -l\-\v -?,\^ . 



(B21) 



As in the /Xo-scheme shown in Appendix IB 1\ the resuh of flB20p can be easily obtained 
from the anisotropic Cgfiv in flS.lOp simply by taking the formal replacements: 



\vi + a,V2 + P, V3 + 7) — > \v + {a + (3 + -f)), 

> 3 ^^'^ p, and Vy-^^ct,V2+/3,V3+j ^ K)+(a+/3+7)5 



/^1,2,3 

where K is the eigenvalues of the volume operator: 

'AirW^^ ( 3^3 A; 



(B22) 



(B23) 



Similarly, the inverse triad operator in the new scheme is given by 



sgn(p) 

¥\ 



\v) 



47r74i/i 



iyli\-v:!_i)\v). 



(B24) 



Classically, C<^ = 87rG'p^/|p|^''^ and this is quantized to the LQC operator C^^\ which when 
acting on the states of 7Y}^°*'^^ yields 



CT\V) 



3 1 

A^) ^B{v)\v)®SnGpl\ 

3^3^ 1 



with 



B{v) 



V + - \v - 1\ 



1/3 



2 /i3 

Consequently, the total Hamiltonian constraint leads to 



(B25) 



(B26) 



"9" 

TlG 
X 



^ + i|i/3_|^_i|i/3r' 

t; + 3|-|t; + l| ?/'(f + 4,0)+ \v - 1\ ~ \v - 3\ V^(t;-4,0) 
t; + 3|-|t; + l| + |f - 1| - |t; - 3| j ^(t;,</ 



where we define 



B{v) ■.= fi^B{v) 



3V3k 



V + - \v - 1\ 



1/3 



(B27) 
(B28) 
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It can be easily inspected that the operator is self-adjoint in L'^iM.Bohr, B{v) dvBohr)- 



Remark: The total Hamiltonian constraint given in ()B27p is almost the same as the master 
equation used in Q]. However, there is a subtle difference. In |8i], B{v) is not introduced and the 
physical Hilbert space is identified as L'^{W^ohry B(v) dvBohr)- Here, on the contrary, the physical 
Hilbert space is L|(RBohr) ) c^wsohr)- 

This difference corresponds to different structures in the WDW theory. In the physi- 
cal Hilbert space in the corresponding WDW theory is L|(R, dv), which is different from 
Ll{R,B{n)d^i). [Note: B{v) = B{fi) := and B{v) B{v) when \v\ > 1.] Here, on 

the other hand, the physical Hilbert space in the corresponding WDW theory is taken to be 
Ll{R,B{v)dv) = Ll{R,B{fi)dn). [Note: B{v) is defined such that B{v)dv = B{fi)dn and we 
have B(v) B_{v) when \v\ ^ 1.] This difference can be understood as two different choices of 
factor ordering both in LQC and WDW theories. Both choices of factor ordering have their own 
advantages. 

The approach used in this paper by introducing B(v) (and B{v) for the anisotropic case) is 
to insist the same WDW theory for both °/i-scheme (or say //o-scheme for the isotropic case) and 
/2-scheme, even though the LQC theory is modified from °/i-scheme to /i-scheme. The advantage of 
this choice is two-fold. Firstly, the expression of the constraint equation is simpler. Secondly, since 
the WDW theory is unchanged, all we have to do in //-scheme is changing variables from to vi. 
Therefore, we can quite straightforwardly "translate" the results we obtained in Section IVII and 
Section [VIII to those in /x-scheme as what we did in Section [VIII CI and Section [VIII Dl However, 
the drawback of this factor ordering is that the gravitational part of the Hamiltonian constraint 
Cgrlv is not self-adjoint in Wgrav^^ 

On the contrary, the choice made in [sl insists Cgrlv to be self-adjoint in Ti^^^ and yields the 
"natural" factor ordering of the WDW theory in accordance with the use of Laplace-Beltrami 
operator in geometrodynamics (see for details). □ 
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